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ABSTRACT 

Cosmological galaxy surveys aim at mapping the largest volumes to test models 
with techniques such as cluster abundance, cosmic shear correlations or baryon acous- 
tic oscillations (BAO), which are designed to be independent of galaxy bias. Here we 
explore an alternative route to constrain cosmology: sampling more moderate volumes 
with the cross-correlation of photometric and spectroscopic surveys. We consider the 
angular galaxy-galaxy autocorrelation in narrow redshift bins and its combination 
with different probes of weak gravitational lensing (WL) and redshift space distor- 
tions (RSD) . Including the cross-correlation of these surveys improves by factors of a 
few the constraints on both the dark energy equation of state w{z) and the cosmic 
growth history, parametrized by 7. The additional information comes from using many 
narrow redshift bins and from measurement of galaxy bias with both WL and RSD, 
breaking degeneracies that are present when using each method separately. We show 
forecasts for a joint w{z) and 7 figure of merit (FoM^^) using linear scales over a deep 
{lAB < 24) photometric survey and a brighter {Iab < 22.5) spectroscopic or very ac- 
curate (0.3%) photometric redshift siirvey. Magnification or shear in the photometric 
sample produce FoM^^ that are of the same order of magnitude of those of RSD or 
BAO over the spectroscopic sample. However, the cross-correlation of these probes 
over the same area yields a FoM^^ that is up to a factor 100 times larger. Magnifica- 
tion alone, without shape measurements, can also be used for these cross-correlations 
and can produce better results than using shear alone. For a spectroscopic follow-up 
survey strategy, measuring the spectra of the foreground lenses to perform this cross- 
correlation provides 5 times better FoM^,^ than targeting the higher redshift tail of 
the galaxy distribution to study BAO over a 2.5 times larger volume. 



1 INTRODUCTION 

Weak gravitational lensing is a unique tool to study the 
large scale dark matter distribution. The correlated shear 
measurements (i.e. shear-shear correlation function) provide 
direct information on the power spectrum of matter fluctu- 
ations (see Bartelmann & Schneider 2000, Refregier 2003, 
Bernstein 2009 for a review), avoiding the problem of bias, 
i.e. of how light traces the mass. The drawbacks are that 
the information is projected on the sky over a very broad 
radial kernel, and that the signal-to-noise is quite poor on 
small scales. For the former, shear-shear tomography can 
be of some help, but the radial precision is still typically 



limited to a few redshift bins (e.g. see Hu 1999, Hu & 
Jain 2004, Refregier et al. 2011). For the latter we would 
need to measure shear for fainter magnitudes and, from the 
ground, this is limited by the atmospheric PSF distortion. 
Cosmic shear also suffers from systematic difficulties, such 
as shape measurement errors, the galaxy intrinsic alignment 
and photo-z errors. Here we will explore the combination 
and cross-correlation of wcak-lcnsing probes with the 3D in- 
formation in the (foreground) large scale galaxy distribution 
in order to complement and extend the constraining power 
of weak lensing and clustering measurements. The galaxy- 
shear cross-correlation can also be used to measure galaxy 
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bias. Once bias is known, we can recover the full 3D matter 
clustering from the 3D galaxy clustering. 

We pay special attention to cosmic magnification as 
an alternative or complementary probe to shear. Lensing 
changes the area of the background image, which can result 
in density fluctuations in the background sources that are 
correlated with the density fluctuations in the foreground 
lenses. The resulting galaxy-galaxy cross-correlation is called 
cosmic magniflcation and contains very similar cosmological 
information to that provided by shear. Magnification will 
be compared with galaxy-shear (e.g. Johnston et al. 2007) 
and shear-shear as a reference (see e.g. Van Waerbeke 2009, 
Bernstein 2009, Van Waerbeke et al. 2010 and references 
therein). Cosmic magnification has long been investigated 
as a potential source for the puzzling galaxy-QSO cross- 
correlation (see Gaztanaga 2003 and references therein). 
Thanks to the improvements in photometric homogeneity 
in the SDSS sample, magnification has been detected in the 
cross-correlation of foreground photometric SDSS galaxies 
with both spectroscopic quasars (Scranton et al. 2005) and 
Lyman break galaxies (Hildebrandt, van Waerbeke & Er- 
ben 2009). In our approach, we envision that, in contrast 
to these current detections, in the near future the lensing 
sample should be the one with the highest radial resolution, 
i.e. spectroscopic or with very good (0.3%) photo-z error. 
For the background sources we do not need such good ra- 
dial resolution to predict the lensing response. This will be 
ideal for 3D power spectrum reconstruction, as will be shown 
here. 

Jain & Taylor (2003) and Bernstein & Jain (2004) pro- 
posed the use of photometric redshift data to produce a tem- 
plate map of foreground galaxies to cross-correlate with the 
induced shear as a function of the background (or source) 
galaxy redshift. They considered the ratios of this cross- 
correlation to recover the geometrical ratios rather than try- 
ing to do a full 3D tomography. Bernstein & Jain (2004) 
called this cross-correlation cosmography. This is related to 
our approach here. The common advantage in both cross- 
correlation techniques is that galaxy-shear cross-correlation 
is less sensitive than shear-shear to systematic errors in the 
PSF correction. The main differences with what we do here 
is that we consider direct cross-correlations, rather than us- 
ing a template, over spectroscopic or very good photomet- 
ric foreground galaxies using narrow bins. We also include 
cosmic magnification, restrict to linear scales (to avoid the 
uncertainties in galaxy power spectrum), and use the full 
signal (rather than just ratios). This allows to recover cos- 
mological information in the 3D power spectrum as well as 
recovery of the geometrical information in lensing efficiency 
(see §3). 

As a new ingredient, we also consider the information 
contained in redshift space distortions (RSD). On linear 
scales, galaxy (peculiar) velocities can be measured by com- 
paring radial to transverse correlations and this provides a 
measurement of the linear growth rate of dark matter fluc- 
tuations. A challenge and opportunity in our approach is the 
need to model galaxy biasing which affects both RSD and 
WL cross-correlations, but in different ways. In real space, 
or for transverse modes, on linear scales, the galaxy auto- 



correlation depends on the square of bias, the galaxy-shear 
cross-correlation only depends linearly on bias and shear- 
shear is independent of bias. This means that their com- 
bination can be used to measure bias as well as the growth 
information (e.g. see Hoekstra et al. 2002 and |2.4| below). In 
redshift space, i.e. for radial modes, bias only affects the den- 
sity growth but not the velocity growth (e.g. see Eq |24| be- 
low), which provides an alternative route to separate growth 
from bias. We will show here that the combination of WL 
probes and RSD results in a very accurate determination 
of bias, also breaking the degeneracy of growth with cosmic 
expansion. Thus, we will combine three different types of 
probes here: 

• Angular clustering from galaxy-galaxy autocorrelation 
in narrow redshift bins (which we call GG) 

• WL from shear-shear (SS), galaxy-shear (GS) or mag- 
nification (MAG, i.e. galaxy-galaxy cross-correlation) 

• RSD, from the ratio of transverse to radial modes. 

The importance of combining redshift space distortions 
and weak lensing information to test cosmological models 
and modified gravity has been highlighted by several authors 
(e.g. Zhang et al. 2007, Guzik, Jain & Takada 2010, Song et 
al. 2011 and references therein). Recently, Reyes et al. (2010) 
used this combination to measure the growth of structure 
from a ratio of galaxy-galaxy auto-correlation with galaxy- 
shear cross-correlation in combination with RSD. Here we 
will see how this can be generalized when we use all the 
information, not only the ratios, and have many redshifts 
bins with much better radial resolution. Potentially RSD 
could benefit from the combination of samples with differ- 
ent biases over the same region of sky (McDonald & Seljak 
2009). This effect is also included in our analysis as we will 
consider bright and faint galaxy populations with different 
biases. But this has little effect on our forecast for the par- 
ticular surveys we are modeling. It might also be possible to 
exploit this further by further division of the galaxy popu- 
lation or direct lensing calibration of biasing (Bernstein & 
Cai 2011). 

Cross-correlation of photometric and spectroscopic 
samples in the same redshift bin can also be used to es- 
timate photometric redshift contamination and determin- 
ing the redshift distribution of the photometric sample (e.g. 
Newman 2008). This technique will be generalized here to 
the cross-correlation between separate bins by inclusion of 
transition probabilities Vij. 

A motivation for our study has been the science case 
for the PAU Survey (www.pausurvey.org) based on the PAU 
Camera (PAUCam, Casas et al. 2010) a photometric camera 
with 40 narrow (lOOA) filters and 6 broad band filters that 
will be commissioned at WHT Telescope in La Palma at the 
end of 2012. PAUCam can map about 2 deg^ of the sky each 
night with all these filters producing an imaging survey to 
a depth of about ias — 24. Because of the narrow band fil- 
ters, galaxies with iAB < 22.5 will have a photo-z accuracy 
of about 0.35% (i.e. about 10 Mpc/h), while the accuracy 
for the rest of the galaxies is about 3% (both with around 
50% completeness for all type of galaxies). The bright sam- 
ple will be much denser ( ~ 15000 galaxies per deg^) than 
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Redshift Space Distortions 
Weak Lensing 
Baryon Acoustic Oscillations 
Galaxy-Galaxy auto-correlation 
Magnification from GG auto and cross-correlation 
Shear-Shear correlations 
Galaxy-Shear correlations 
combination of SS, GS, GG and MAG 
Fain sample of galaxies 22.5 < iab < 24 
Bright sample of galaxies iAB < 22.5 
Combination of independent F and B samples 
Cross-correlation of F and B over same area 



Table 1. Notation and acronyms used in this paper. 



any spectroscopic surveys to the same depth. Unhke most 
spectroscopic surveys, which target predetermined galaxies, 
PAU will contain all the objects in the surveyed area to 
a given brightness. In our initial studies for PAUCam, we 
only explored BAO with iAB < 22.5 LRG galaxies (Benitez 
et al. 2009). However, as we discuss next, this is a rather 
restrictive use for the PAUCam data. A similar approach 
has been taken by several spectroscopic surveys which plan 
to measure BAO at high redshifts by selecting appropriate 
spectroscopic targets out of a given parent photometric cat- 
alog. Here we will consider forecasts using both the faint 
photometric and bright (quasi) spectroscopic samples and 
we will show that one can get significantly better cosmologi- 
cal constraints when using probes based on the combination 
of WL and RSD. 

We will consider both a modest 200 deg^ PAU-like Sur- 
vey and a more ambitious 5000 deg'^ survey. Such a survey 
could also result from the Dark Energy Survey (DES), which 
will image 5000 deg'^ in five passbands to 24th mag over 5 
years starting in 2012, producing weak lensing shape mea- 
surements for about 200 million galaxies. A massive follow- 
up spectroscopic survey over a substantial part of the DES 
footprint, e.g., using multi-fiber spectrographs such as Big- 
BOSS or DESpec, would enable RSD measurements of the 
lensing population and the implementation of this cross- 
correlation technique. Unlike BAO or supernovae, WL and 
RSD can also provide very valuable information on cosmic 
growth history and this is a key ingredient to understand 
the physics of the accelerating Universe. 

This paper is organized as follows. In §2 we present the 
modeling, methodology and different assumptions used. We 
include a subsection on how we model galaxy bias and an- 
other one on modeling of redshift space distortions. Section 
3 describes the approximations we make for weak gravita- 
tional lensing. The reader familiar with these techniques can 
directly jump to §4, where we present our fiducial surveys. 
In §5 we present results of our forecast as a function of the 
different ingredients , to disentangle the different contribu- 
tions to the FoM. This section could also be skipped if the 
reader is not interested in such details. In §6 we present the 
main result in this paper, i.e. the comparison of forecasts 
for different surveys. We finish in §7 with some conclusions 
and summary of the main results. Tabl^ summarizes the 
notation that will be used in this paper. 



2 MODELING 

In this section we will introduce the different assumptions 
and modeling used in this paper. Section 2.1 and 2.2 intro- 
duce the cosmological model and the parameters that we 
want to study. In Section 2.3 we present the Fisher Matrix 
approach and the figures of merit that will be used to com- 
pare experiments. Sections 2.4 and 2.5 present and justify 
the models for galaxy bias and redshift space distortions. 



2.1 Growth and Cosmic History 

The cosmic expansion history, a — a(t) or H = H(t), in a flat 
Friedman Lemaitre Robertson Walker (FLRW) background 
with matter density pm and dark energy (DE) equation of 
state w — Pde / pDE, can be written as: 



SttG 



{pm + Pde) 



(1) 



Q.DEa 



-3(l + w) 



where Q,k = 1 — Q.m — ^DE measures deviations from flat cur- 
vature k = and we have neglected radiation. Our fiducial 
model corresponds to ACDM: a fiat universe with w = —1, 
so that the DE density is constant with redshift z. We will 
explore how well our different observational probes can con- 
straint w and its variation w = 1(7(2:) Parameters wq and 
Wa are used to characterize the evolution of DE equation of 
state (Chevallier & Polarski 2001, Linder 2003): 



w{z) = Wo + Wa{l — a) = Wo + WaZj (1 + z) 



(2) 



According to General Relativity (GR), given this cosmic 
history, the equations that determine the growth history, i.e. 
the cosmic evolution of the linear density contrast 5, are of 
the form (Peebles 1980; Bernardeau et al. 2002) 



with the solution 



& + 2HI, = 47rGp„5 



D(a)(5(0) 



(3) 



(4) 



where the growth factor D{a) depends on the expansion his- 
tory Ilia) (through wia)) and on S7m(a). Any discrepancy 
found between the observed growth and the growth D pre- 
dicted for a given expansion history H can be use as a test 
for modifications to GR or variations on the cosmological 
model. This linear growth can also be characterized by its 
derivative, the velocity growth factor: 



dlnD 
d In a 



(5) 



where 7 is the gravitational growth index (see Linder 2005). 
So when normalized to _D = 1 today, then 



D(a) 



■ exp 



din a /(a) 



(6) 



^ When w = w{z) we need replace w in floEO- in Eq[l] 

by the corresponding integral over redshift. 
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For GR with DE equation of state w, Linder (2005) finds 
that to a good approximation: 

3{w - 1) 



7 ■ 



Gw — 5 



(7) 



This reduces to the well known result 7 ~ 0.55 for A.CDM 
(w — —1). Thus for DE models in GR, measurements of f{z) 
can be used to estimate ^(2), independently from measure- 
ments of H{z). Other cosmological models have a different 
relation between D{z) and H{z) (e.g. see Gaztanaga & Lobo 
2001) which results in different effective values for 7. For 
example, in the DGP model 7 ~ 0.68 (Lue, Scoccimarro 
& Starkman 2004). More generically, an independent mea- 
surement of 7 and w can be related to a measurement of 
time variations of Newton's constant G (e.g. see Eq.32 in 
Pogosian et al. 2011), which provides a direct test of GR. 
This will be our approach here: we will check how well we can 
measure separately w and 7 with different probes, assuming 
they are independent of each other. We will not make use of 
GR to relate cosmic growth to cosmic history (i.e. we will 
not use Eq[7] or Eq|3|, but rather we will use a paramet- 
ric description in terms of independent parameters: wq, Wa 
and 7 (plus the standard cosmological parameters in Eq[9] 
below). This can be used to separate GR from other theo- 
ries of gravity and learn about the nature of the accelerated 
expansion. 



Note that for completeness we include as = cr8(0), the 
global normalization of the amplitude of density fluctuations 
at z = 0. Our goal is not to measure erg but rather the evo- 
lution of this amplitude with time, i.e. D{z), which is given 
by 7. In our analysis, as is dominated by the priors above 
(around 1%) as our probes are more sensitive to relative 
than to absolute variations. One could also try to measure 
as and other parameters, such as neutrino masses, by in- 
cluding the amplitude of CMB temperature fluctuations at 
z=1100 rather than using priors of as (i.e. at z — 0). This is 
a very different approach and involves other complications. 
In such case 7 and as could be strongly correlated. But this 
is not our approach here. Our goal is to test the best way 
to measure cosmic and growth evolution starting from some 
priors at z — 0. 



2.3 Fisher Matrix and Figures of Merit 

Throughout this paper we will employ the formalism based 
on the Fisher Matrix (FM, e.g. Tegmark et al. 1998): 



2.2 Fiducial Cosmological model 

To characterize cosmic history and the linear matter power 
spectrum, we use the standard set of 8 cosmological param- 
eters (e.g. see Komatsu et al. 2011): 



WOjWa, h, Us, 0,m, f^B, ^DE, as 



(8) 



In addition, the growth index in Eqj5]is allowed to vary in- 
dependently of w(z) to characterize the growth history and 
to test Modified Gravity separately from DE evolution. To 
model bias we include an additional set of independent bias 
parameters as described later on in j ]2.4| Thus we have a 
total of 9 cosmological parameters together with some nui- 
sance bias parameters 6^: 



(wo, Wa,h, ns,ilmMB,^lDE,as, hi, 7) . 



(9) 



Note that in general fifc = 1 — Q.m ~ Q,de differs from 
zero. For the fiducial values we use wo — —l,Wa — 0,h — 
0.7, Us = 0.95, fi™ = 0.25, fis = 0.044, f^ois = 0.75,^8 = 
0.8, 7 = 0.55 (corresponding to its value in GR). Unless 
stated otherwise, results are always presented with priors 
from Planck and Stage-II Supernovae (SN-II)j^ but with no 
priors in 7 or bias bi. 

2 We use the Dark Energy Task Force (DETF) 
Planck and SN-II priors Fisher Matrices given in 
|http:/ /ww w. physics. ucdavis.edu/DETFast/ with file names 
planckfish and SN-II. The Planck prior on trg assumes a GR 
model for linear growth but there are several other ways to 
measure as that could produce similar priors without such 
assumption. The actual priors and how they were obtained 
are somewhat irrelevant here, they just represent a baseline to 
compare different experiments. 



J2 E ^ 



£ or fc ij,mn 



dPti 



n 



(10) 



where are the cosmological parameters, nuisance param- 
eters or additional quantities we want to forecast and dj 
are the observables with Q covariance. The indices ij refer 
here to each pair of redshift bins and i or k refers to angu- 
lar scales: multipole or Fourier modes. The elements of F^J 
are the covariance of and p^. For a given experiment, the 
FM approach is a mapping (or transformation) of the covari- 
ance in the observables to the covariance in the quantities 
we want to measure. 

Here we will model the weak lensing (WL) probes, i.e. 
shear-shear, galaxy-shear and galaxy-galaxy angular corre- 
lations, in real space and using the Limber approximation 



(see i 3.2 1. In this case the observables Cij will be the power 
spectrum of the angular cross-correlations between two red- 
shift bins. Redshift bins will be taken to be independent, 
as they are separated and non overlapping in space. This 
means that there are no intrinsic radial correlations and we 
only consider transverse modes for angular clustering and 
WL probes. For redshift space distortions (RSD) we will 
use the ratio of the amplitudes of the 3D power spectrum 
as we change from transverse to radial modes within a sin- 
gle redshift bin, i.e. fj, in Eq |24[ We neglect the covariance 
between these ratios and the transverse angular power spec- 
trum. When RSD is combined with BAO we also neglect 
the covariance between the RSD ratios and the BAO wig- 
gles. Hence the combination of RSD with the WL probes or 
BAO is given by the addition of the corresponding Fisher 
matrices. 
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2.3.1 Non-linear scales 

The FM depend on scales |^ as the summation is over all 
scales that contribute to the analysis in Eq. (10 1. If we fix 
a minimum scale Rmin defined by the validity of "linear" 
theory today, the redshift evolution of Rmin can be obtained 
through this relation: 



0-{Rmin, Z) 



(11) 



where (j{R,z) = / Pi,(k, z)W'^ {kR)cfk is the amplitude of 
fluctuations on scale R, W{x) — exp(— a;'^/2) is a Gaussian 
smoothing kernel and Pi,{k,z) is the linear matter power 
spectrum. These are the same criteria adopted by White et 
al. 2009. The maximum fe-mode is then given by 



, Rmin\0) ,7.^ — 1 
kmax MpC 



(12) 



Rmin (■2') 

with kmax = 0.1 to "normalize" this scale to 0.1/iMpc~^ at 
z = 0. The corresponding maximum £-mode is given by 



+ 1/2 = kn 



)r{zi) 



(13) 



where r is the transverse comoving distance to the redshift 
bin i under consideration. In the forecast, we consider cross- 
correlations to scales up to the Imax of the closest bin (i.e. 
the smallest Imax of the two bins). For the maximum scale 
we use kmin = 0.001 ^Mpc~^ for all z and Imin = 50 to 
avoid inaccuracies in the Limber approximation!^ 

These conditions, in Eq.( 11 ]-( 12 1, yield very similar val- 
ues for kmax to those found by studying the limitations of 
perturbation theory by Jeong & Komatsu 2006 (see also 
Jeong & Komatsu 2009). In this analysis we will only use 
linear theory for P{k) and Gaussian initial conditions. 



2.3.2 Figures of Merit 

In the FM approach the marginalized 1-sigma variance in 
parameter pv is given by the diagonal element of the inverse 
of the FM: a^{p^) — [F~'']^^. When optimizing an experi- 
ment with more than one parameter, you need to take into 
account their covariance. We will then cast our results in 
terms of: 



FoMs 



det[F- 



(14) 



where S is the sub-space of parameters we are interested 
on. For just one parameter, this is the inverse error, for two 
this is proportional to the inverse area included within the 
1-sigma error ellipse (Albrecht et al. 2006), for three this 
is the inverse volume within the 1-sigma error ellipsoid, and 
so on. In this paper, we will focus on three different figures 
of merit: 

• FoM„ : Constraints on {w(),Wa) marginalizing over the 
remaining cosmological and bias parameters, similar to the 

^ We employ spherical coordinates in the formalism of weak lens- 
ing (i.e. £-modes) and f modes for RSD (i.e. fe-modes). 
* Notice that error-bars are larger at large scales, hence this 
choice has little impact on FM constraints. 



DETF FoM (Albrecht et al. 2006). But note how here we 
allow 7 to vary, while in Albrecht et al. (2006) this figure of 
merit is quoted for a fixed value of 7 ~ 0.55 corresponding 
to GR. 

• FoM.y : the inverse of the achievable 1-sigma error in 7 
after marginalization over the remaining parameter space. 

• FoMu,-, : Joint constraint on (wo,uia,7), a natural ex- 
tension to the dark energy FoM of the DETF. 

We find that there is little covariance between 7 and 
(woyWa), which means that FoMu,-y ~ FoM^x FoM^,. Note 
that the above definitions have different dimensions since 
they involve different numbers of parameters. To have a di- 
rect comparison of the (geometrical) mean error in the dif- 

9 2/3 

ferent cases we should compare FoMu, to FoM.^ and FoM^,-, . 



2.4 Galaxy Bias and its cross-correlation 

In this section we will justify how we model galaxy bias. This 
is a key ingredient in our approach as both RSD and WL 
cross-correlations depend on bias. We will show that we do 
not need to include a correlation coefficient between galaxy 
and matter cross-correlations, and argue that we only need 
to allow bias to vary on time scales larger than Aa ~ 0.1. 



2.4.1 Bias stochasticity and scale dependence 

In the local bias model (Fry & Gaztanaga 1993), the 
smoothed galaxy fiuctuation 5g at a point a; is a function 
of the matter field 5m at the same point. For small fiuctua- 
tions we can Taylor-expand this function and approximate 
the local bias model as: 



3g{x, z) = b{z) 5m{x, z) 



(15) 



The galaxy-galaxy (auto) correlation £^gg =< Sg5g > and 
galaxy-mass (cross) correlation =< 5g8m > are then: 

^gg(s,z) = fe^(z) Cmm(s, 2) OC h^{z)D'^{z) 
^gm(s,z) = b{z)^mmis,z) OC b{z) (z) (16) 

where (,mm{r) is the matter-matter correlation and s = 
|a;2 — is the separation between pairs and the second step 
uses the linear gravitational growth. Under these linear and 
local assumptions one can combine ^gg and (^gm to break the 
degeneracy between b and D. A more general parametriza- 
tion of bias is based on the cross-correlation coefficient, r, 
defined as: 



(17) 



If instead of using the local bias model we assume an effec- 
tive bias P = ^gg/^mm that is constant as a function of scale 
we find: 



\/ ^ga / ^rnm — 
^gm/^mm — 

igg/^gin ~ 



b(z) 

b(z) r{z) 

M 

r{z) 



(18) 
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where r{z) can difTer from unity (e.g. see Tegmark & Peebles 
1998, Pen 1998, Dekel & Lahav 1999, Casas-Miranda etal 
2002, Seljak & Warren 2004, Bonoli & Pen 2009). In the 
local model r{z) = 1 and all these ratios are just equal. It 
is still possible to measure D{z), b{z) and r{z) using more 
observables, such as ^mm from shear-shear or P{k, ^) from 
redshift space distortions. 

We will show next that on large scales, when b is con- 
stant as a function of scale, r has to be close to unity and 
there is no need to introduce this parameter. Let us add non- 
linearities and a stochastic component e{x, z) which could be 
caused by shot-noise or non-local perturbations: 

6g{x,z) = bi{z) S,n{x,z) + b2{z) Sl^(x,z) +e{x,z) (19) 

By definition, e is not correlated to (5,„. Otherwise, the cor- 
related part will just renormalize bi and 62. Then 



= 61 



+ + 2bib2^3 + blU 



(20) 

> and 



where we have introduced =< ee >, ^3 =< (5„ 
^4 =< (5m 5m > for the 3 and 4 point correlations. All these 
new correlations have a different scale dependence on the 
separation s than ^gg(s). For Gaussian initial conditions: 
^3 ~ and ^4 ~ 5mm (see Bernardeau et al. 2002). On 



large scales, s > lOMpc/h, where 5n 



< 1: 



C99 

5gm 



t'l smm mm 



bU. 



bl 5r7 



\/l + UbUr, 



26?5" 



(21) 
(22) 



In general (,e{s) « 5mm (s) and the above equations repro- 
duce the local bias results in Eq |16| with b — bi, even when 
the stochastic e can be shown to be quite large, even larger 
that Sm (see Fig.l in Manera & Gaztanaga 2011). When 
5e(s)/5mm(s) is uot negligible then we have that r is scale 
dependent, since both 5e and 5mm depend on the pair sep- 
aration s. But, in such a case the effective bias, defined as 
1? = ^gg/^mm, will also depend on scale. In other words, 



whenever the effective bias b is independent of scale, Eq 21 
indicates that we can neglect the sc ale d ependent terms, i.e. 
5E/&^5mm and then, according to Eq 22 we force r ~ 1. This 
result is well reproduced in halo bias simulations on scales 
s > 20Mpc/h (see Fig. 4 and Fig. 13 in Manera & Gaztanaga 
2011 and next sub-section). 

Seljak & Warren 2004 and Bonoli & Pen 2009 found 
that r tends to unity on large scales in halos in numerical 
simulations but deviations from r = 1 are significant even on 
the largest scales. It should be noted that these references 
studied r{k) in Fourier space for halos and not in configura- 
tion space. In Fourier space r is also subject to shot-noise, 
because r is build from the power spectra < 5{k)^ >, while 
shot-noise cancels for the correlation < 5{r)5{r') >. In our 
forecast we use fourier and harmonic space, but we want 
to model galaxies and not halos. For galaxies or optimaly 
weighted halos shot-noise is less important and r is closer 
to unity (see Scherrer & Weinberg 1998, Matsubara 1999, 
Hamaus etal 2010, Cai & Bernstein 2011). Correcting for 
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Figure 1. Bias in mock MICE galaxy simulations. Mean and 
(2-sigma) errors correspond to a 5000 deg'^ survey. Top panel 
shows the evolution of the effective bias defined in different ways: 
y/igg/imm (squares), igg/S,gm (circles) or S,gm/S,mm (triangles). 
The agreement is excellent: we include a small horizontal dis- 
placement to be able to see the different symbols. Continuous (or 
dashed) line shows the best fit bias evolution model: b(z) based on 
linear inter/extrapolation between 4 (or 2) free bias parameters 
b(zi) spaced at regular scale factor values a = 0.8, 0.7, 0.6, 0.5 (or 
a = 0.7, 0.5). Bottom panel shows the cross-correlation coefficient 



shot-noise of halos is possible but complicated as halos have 
a exclusion region which results in a sub-Poisson correction 
(see Casas-Miranda etal 2002, Manera & Gaztanaga 2011). 
Below we will show that indeed r is very close to unity in 
numerical simulations for all redshift. In §5.2.3| we will test 
how a stochastic bias impact our result. 



2.4-2 Bias evolution 

Figure [l] illustrates the above modeling as a function of 
redshift. Galaxy mocks are built from halos in the light- 
cone of a MICE simulation of side L — 3072Alpc/h and 
A*' = 2048^ particles (see Fosalba, Gaztanaga, Castander & 
Manera 2008, Crocce, Fosalba, Castander & Gaztanaga 2011 
for more details on the simulation). The large size of this 
MICE run allows to build a light-cone sampling one octant 
of sky without repetition to z = 1.5. Halos are identified and 
weighted with a simple HOD (Halo Occupation Distribu- 
tion) prescription. The number of galaxies in a halo of mass 
M > Mruin is l + (2oi?£— )° ^ith Mmin ^ 10^^M^„„//i aud 
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Q ~ 1 (e.g. see Scoccimarro et al. 2001) Evolution with 
redshift is all given by the halo mass evolution. We divide 
the light-cone simulation in redshift bins of width Az ~ 0.1 
and estimate the correlations ^mm, £,gm and ^gg in each bin 
as a function of separation s. Three different bias parameters 



are then fitted to the ratios in Eq 18 on scales s > 30Mpc/h 
where we find that these ratios are constant within the er- 
rors. Results for the different ratios agree well as shown in 
the top panel of Figjl] These results indicates that r ~ 1, 
as confirmed by the bottom panel. Errors are from a 100 
jackknife angular patches which we have checked agree well 
with subsample errors (e.g. see Norberg et al. 2009 and ref- 
erences therein). The error-bars are similar for the difi'erent 
cases because sampling variance cancels to a good extend in 
the ratios. The bottom panel shows that is indeed very 
close to unity in agreement with the arguments above. The 
mean rms error in r is only 0.2% and the rms deviation from 
unity is less than —0.1%. Similar results are also found for 
different choices of HOD weights. 

Note that because of the large volume and density of our 
simulation, the evolution of bias is known here to better than 
1%. Halo model predictions do not yet reach these levels of 
accuracy (Manera, Shech & Scoccimarro 2010; Manera & 
Gaztafiaga 2011), but could potentially be calibrated with 
simulations to achieve similar accuracies. As shown in the 
figure, the values at different redshift are strongly correlated, 
indicating that only a few parameters are needed to char- 
acterize the evolution of bias. The characteristic time scales 
for bias evolution in our samples seem to be Aa > 0.1, corre- 
sponding to t > IGyr, which is typical of galaxy evolution. 
This is also in agreement with recent bias measurements 
(e.g. Coupon etal 2011), which find a smooth dependence of 
observed bias with redshift for different galaxy luminosities 
and color, in good agreement with HOD. These results cor- 
responds to iAB < 22.5, similar to the sample that will be 
studied here. Systematics and selection effects do not seem 
to introduce higher frequency variations in the recovered bi- 
ases. 

Following these indications, we model b{z) by including 
free bias parameters bi at some fixed redshifts, with lin- 
ear interpolation to other redshifts. We choose a constant 
spacing in scale factor Aui = 0.1 to fix the interpolation 
positions in b{z). Given the redshift distribution of galax- 
ies in our sample we choose 4 points as the interpolation 
locations: fli = 0.8,0,7,0.6,0.5, corresponding to redshift 
Zi — 0.25,0.43,0.67,1.0. Fig[l] shows, as continuous line, a 
fit to such bias model. The fit has a per degree of freedom 
of 0.26 which is "too good". This indicates that this model 
has too many free parameters. If we use 2 bias parameters, 
instead of 4, at — 0.7, 0.5, we find a per degree of free- 
dom of 1.2 (dashed line in the Figure). This indicates that 
2 parameters is probably not enough, given the small error- 
bars. So 3 parameters seems to be the right number for large 
samples (5000 deg^) and 2 parameters could be enough for 



^ This assignment can also introduce additional stochasticity but 
its contribution to vanishes if this stochasticity is not spatially 
correlated between halos 



smaller area surveys (ie 200 deg^) which have larger errors. 
Here we will always use 4 parameters to be on the safe side. 

In summary, we will use 4 biasing parameters with no 
priors as default, but also show how results change as a 
function of priors. Both from HOD modeling and from its 
comparison with observations we already have strong priors 
on what is the variation of bias for a given galaxy sample. 
These priors can also be estimated from the same data we 
are considering here. Recall that we will only be using linear 
scales in our forecast, while HOD modeling can take advan- 
tage of the data on smaller non-linear scales to constraint 
halo mass and HOD. Moreover, biasing could also be mea- 
sured and constrained with other techniques, such as higher 
order correlations in the same galaxy sample (e.g. see Gaz- 
tanaga, Norberg, Baugh & Croton 2005, Sefusatti, Crocce, 
Pueblas & Scoccimarro 2006 and references therein). 



2.5 Redshift Space distortions (RSD) 

The measured redshift distance s to a galaxy differs from 
the cosmological distance r by its peculiar velocity v along 
the line-of-sight. These displacements lead to redshift dis- 
tortions. On large (linear) scales the dominant effect is due 
to coherent bulk motions induced by gravity. Mass conser- 
vation implies that the velocity divergence is 6 = V.w = ~5, 
which on linear scales results in 6 — —fS (see Eqj5|. Thus 
in the line-of-sight direction fiuctuations are distorted by a 
factor (1 -I- /) resulting in a squashing effect in the 2-point 
redshift correlation function (Kaiser 1987). At small scales, 
random velocities inside clusters of galaxies produce a radial 
stretching pointing towards the observer, known as fingers 
of God (FOG). This have very little effect on linear scales 
considered here (Hikage, Takada & Spergel 2011). 

The top left panel of Fig[2] shows the squashing effect 
in the amplitude of the 2-point correlation as a function 
of radial (vertical axes) and perpendicular (horizontal axis) 
separation (Gaztafiaga, Cabre & Hui 2009). The squash- 
ing effect is quite clear in the inner regions of 20-50 Mpc/h 
and it also produces a large region with negative correlation 
between 50 and 100 Mpc/h which is characteristic of RSD 
(it is not present in the transverse direction). This effect al- 
low us to measure / (Cabre & Gaztafiaga 2009). One can 
also use this information to measure H{z) from the BAO 
ring position at ~ 100 Mpc/h (see Gaztanaga, Cabre & Hui 
2009, Matsubara 2004) and put constraints on DE equation 
of state (Gaztanaga, Miquel & Sanchez 2009). 

The remaining panels display the same measurements 
done over a distribution of objects with photometric red- 
shifts of increasing photo-z error (as labeled) . The distortion 
pattern due to coherent infall for the case of = 0.003(1 -I- 
z) (top right panel) is remarkably similar to the one in the 
original redshift sample. And it is precisely this anisotropic 
signal that we will employ to disentangle bias and growth of 
perturbations. For most purposes ~ 0.003(l-|-z) is almost 
equivalent to having a (spectroscopic) redshift sample. This 
will be quantified better later. Notice how this agreement 
degrades rather quickly with az- 

In the large-scale linear regime and in the plane-parallel 
approximation (where galaxies are taken to be sufficiently 
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Figure 2. Top left panel shows the ^(Tr, cr) correlation in the Kaiser model with no photo-z error, i.e. cr^ = (from Gaztafiaga, Cabre & 
Hui 2009). The correlation is clearly squashed in the radial direction with a region of negative correlation (in blue) between tt = 50— 100 
Mpc/h. Top right panel shows the same model but with a photo-z degradation of CTz = 0.003(1 + z), corresponding to the PAU Survey. 
The difference is small and is mostly confined to small radial scales. Bottom panels show how the results are degraded as we increase CTz 
to 0.006 (left), 0.016 (center) and 0.033 (right panel). As the photo-z increases the radial squashing disappears, turning instead into a 
radial elongation. Note also how the region of negative correlation vanishes as we increase the photo-z error. 



faraway from the observer that the displacements induced 
by peculiar velocities are effectively parallel), the distortion 
caused by coherent infall velocities takes a particularly sim- 
ple form in Fourier space (Kaiser 1987): 



5.(fc,M) = (l + .//i')5(fc) 



(23) 



where /i is the cosine of the angle between k and the line- 
of-sight, the subscript s indicates redshift space, and f{z) 
is given by Eq. ([5|. If we assume that galaxy fluctuations 
are linearly biased by a factor b relative to the underlying 
matter density 5 (i.e. Sg = 6 5) but velocities are unbiased, 
then 



(24) 



where 5g are the measured galaxy fluctuations in redshift 
space. However in practice what we measure is the rms value 
a{h+ IJ? f), where cr is the r.m.s. amplitude of fluctuations at 
some scale. Hence by fitting the ^ dependence we can have 
independent measurements of 6 cr and f a &s & function of 
redshift (or the ratio /3 = f /h, but we can not separate / 
from both cr or 6). This means that f{z) can not be measured 
unless we flx the normalization of 5{z). Combining RSD with 



weak lensing (see below) or higher order correlations in the 
same galaxy sample (e.g. see Gaztanaga, Norberg, Baugh 
& Croton 2005, Sefusatti, Crocce, Pueblas & Scoccimarro 
2006 and references therein) we can break this degeneracy. 
Nonetheless having an estimate for f{z)a(z) (independent 
of bias) can be as valuable as just having /, as both are ways 
to constrain D(z) and therefore 7. 

To implement RSD constraints on the growth rate / and 
bias we follow the approach of White et al. (2009). First we 
write the cross power spectrum of galaxy samples A and B 
at the same redshift bin i as, 

P^Aik, = {bA + fti'){bB + //i') G., (fc, ^i) PUk) (25) 

where A and B are indices for different galaxy types and Gz 

G.(fc, fi) = exp [-{l/2)k^fi\a\ + al)] (26) 

accounts for the radial smearing due to photometric redshift 
uncertainties (and possibly the small scale peculiar velocity 
dispersion). The factor a relates to the photo-z error by 
cr = azc/H{z). In Eq. (25 1 the bias b, the growth rate / and 
the photometric uncertainty cr are functions of redshift Zi. In 
turn, the linear power spectrum Pl is determined up to an 
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overall normalization, say crs{z) = D{z)asiO), leading to the 
well known degeneracy between f{z) or b{z) and aaiz). The 



FM corresponding to a given redshift bin i is (see Eq. (lOl), 



E 



Vo (fk I dP_ 

(2^)3 I dq^ 



C^V ( ^ 1 (27) 



where Vo is the total volume of the redshift bin and X,Y = 
(ffiPij 5152, (?2<72) in the case of two galaxy types or just X — 
gl for only one type. The derivatives in Eq. ( 27 1 for just one 
population are simply given by, 



dP _ 

db ^ b + 



^ P{k) and 



df b + f^i- 



■P{k), 



while the covariance matrix is, 

Cov[P(fe), P(fc)] = 2P^{k)N^{k) (28) 

with N{k) = 1 + (nP(fc))^^ and n the galaxy number den- 
sity. It is in this error term that the damping due to pho- 
tometric uncertainty has an impact by increasing the shot- 
noise contribution to N{k) and therefore the total error in 
P(k). Expressions for the derivatives and covariance of two 
or more galaxy types can be readily obtained (see White et 
al. 2009). 

The combination of parameters that are constrained 
with RSD are = b{z)Gs{G)D{z) and f{z)(Js.{Q)D{z). In 
turn these depend on the set of cosmological parameters 
defined in Eq. ^ through Eqs. ( 1|5|6 1. To be able to join the 
constraints from RSD and WL we transform the variables 
by: 



(29) 



where M = dp/dq. 



As explained in S 2.2 the value of cjs{Q) is dominated by 



priors so that we only measure relative amplitudes here. As 
we do not measure the shape of -P(fc) either, in practice our 
forecast is equivalent to a measurement proportional to the 
ratio of the radial and transverse directions. 

When we combine RSD with WL measurements we as- 
sume that these measurements are not correlated. This is not 
quite correct when the measurements are over the same area 
(FxB case below). To avoid the complication of estimating 
the full covariance of these two observables (i.e. P{k) with 
d) we will assume here that they are 100% correlated and 
make sure that we only count once the same modes (i.e. we 
use only independent measurements). This is a conservative 
approach as the true covariance will always be smaller. The 
total number (A'^r) of fc-modes in the RSD case is: 

3 



dNrik) 



2-Kk dkdjj, 



(30) 



where is the volume of the survey region (e.g. within a 
redshift bin). This can be written as the product of purely 
transverse Nt and purely radial Nr modes: 



dNrik) = 
dNr{k) = 

dNt{k) = 



dNt{k) dNr{k) 
L 

2%, 



(31) 



kdfi 



2tt 



2-Kkdk {pL = 0) 



In WL (or angular clustering) we will use the Ci an- 
gular power spectrum which only measure purely transverse 
modes. This is easy to verify. If we use k = l/x (where x 
is radial distance and I is the angular multipole) we find 
that Nt{k) ~ N{t) ~ U, as in the Ci calculation in Eq|6Tl. 



So we can avoid double counting by subtracting the purely 
transverse modes from the RSD calculation: 



Nrik) = Nrik) - Nt(k) = Nrik) 1 



Nr{k) 



(32) 



where Ay is the number of independent modes that can be 
used for RSD forecast when combined with WL (ie for the 
FxB case). This correction is small as Nt{k) » 1. In our 
case this can only affect the FxB combination and we have 
checked that this correction is always smaller than 30% in 
the final FoM. 



3 WEAK LENSING 

Weak lensing produces both a coherent distortion of galaxy 
shapes (known as shear distortion) and coherent magnifica- 
tion (MAG) of area that can be a source of density fluctu- 
ations. The latter is known as magnification bias or cosmic 
magnification (see e.g. Gunn 1967, Narayan 1989, Broad- 
hurst et al. 1996 and references therein) which introduces 
both a correction to the angular galaxy auto correlation 
function (Villumsen, Freudling & da Costa 1997; Loverde, 
Hui & Gaztanaga 2008) and also to the cross-correlation of 
galaxies at different (including disjoint) redshifts (Moessner 
& Jain 1998; Hildebrandt et al. 2009). Magnification bias 
also changes the shape of the angular auto correlation func- 
tion in RSD (Matsubara 2004; Hui, Gaztanaga, & Loverde 
2007; Hui, Gaztanaga, & Loverde 2008), but these are rela- 
tively small effects, so we neglect then here from now on. 

Measurements of cosmic magnification can provide com- 
parable S/N to that of cosmic shear without the need to 
measure galaxy shapes (see e.g. Van Waerbeke 2009 and 
references therein). With spectroscopic or good photomet- 
ric redshifts, such as those in PAU, one can further explore 
3D lensing tomography by means of the angular cross cor- 
relation between galaxies in different narrow redshift bins, 
as will be explored here. 

We follow the approach of Hu & Jain (2004 and refer- 
ences therein) who presented forecasts for shear and galaxy 
cross-correlations. As new ingredients, we will also consider 
the case for weak lensing MAG and the use of narrower red- 
shift bins. 

As we will show below, i.e. Eq |56| and Eq |59| both 
galaxy-galaxy and galaxy-shear cross-correlations of narrow 
redshift bins produce a direct measurement of the power 
spectrum of the foreground (lensing) distribution. This al- 
lows for a 3D reconstruction of the power spectrum, with 
a radial resolution that is much better than possible with 
the shear-shear reconstruction, i.e. Eq |58[ The smaller the 
bin width, the better the reconstruction. In our forecast 
we are limited radially by the accuracy of the Limber ap- 
proximation, which breaks for very narrow bins. To simply 
the analysis we also want to be able to neglect the intrin- 
sic cross-correlations of separate bins. To achieve this, we 
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find that we should use redshift bins which are larger than 
about 0.014(l+z) (i.e. > 60 Mpc/h) for which these approx- 
imations are adequate. Better results could in principle be 
obtained by using smaller bin widths, but this is left for a 
future analysis. In this sense our results are conservative. 



3.1 Magnification & Sliear 



Magnification jj, is defined as: 



M = 



1 



detA (1 - k)2 - |7|2 



(33) 



where A is the Jacobian matrix for the lensing transforma- 
tion (see e.g. Bartelmann & Schneider 2000). In the weak 
lensing limit, fluctuations in magnification fi, convergence 
K and shear 7 are closely related. Fluctuations in (E-field) 
shear e and convergence are equal, 5^ = Se (see e.g. Hu & 



Jain 2004) and, according to Eq 33 they are half as large 
as magnification Sfj, — 25k- As k can be obtained from shear 
measurements, we will also use k to refer directly to shear. 
But we should bear in mind that some shear estimators are 
given in terms of Sn/2 (e.g. the iCosmos software of Re- 
fregier et al. 2011) and other combinations. In our FM ap- 
proach this will only be relevant when introducing the scale 
for the intrinsic noise ellipticity. 

We will focus here in convergence, 5k, as reconstructed 
from galaxy shapes measurements (that we will call shear), 
and magnification, 5^, as estimated from fluctuations in 
galaxy number density counts. Magnification changes the 
area of the background sources behind lenses, this induces 
a background fluctuation 5g ~ —5^ which is correlated with 
the foreground galaxy population. Additionally, background 
magnitudes are also aflected inducing additional galaxy den- 
sity fluctuations in Sg across the sample magnitude limit. 
Adding both contributions gives 



5g = (2.5s - 1)5„ ~ (5s - 2)5k 



(34) 



where s here is the slope of the galaxy number counts at the 
flux limit. Dust extinction in the lenses can also produce 
signiflcant fluctuations. Menard et al. (2010) have shown 
that at sufficient large wavelengths, i.e. / and Z bands, dust 
extinction becomes negligible and the change in magnitude 
is dominated by magnification. In our analysis we will use 
I magnitudes to select galaxies and neglect dust extinction, 
although this might not always be a good approximation 
(see Fang et al. 2011). 

Weak lensing convergence in bin j is given by the pro- 
jected matter density (J^. in all foreground redshifts i < j: 



(35) 



where pij < 1 is a geometrical WL weight that will be in- 
troduced later on and 9 are sky positions. We then have 
the following relation between observed galaxy fluctuation 
at background bin j, i.e. 5g. , and the foreground DM distri- 
bution: 



bj5„ 



,(e) + e,(e) 



^Pij5rr 
i<3 



m (36) 



where pij — (5s — 2)pij and we have used the linear lo- 
cal bias with stochasticity e, deflned in Eq |19[ If we neglect 
e we can use the observed 5g{9) and Sk{0) maps and the 



above Eq 36 and Eq 35 to reconstruct both 5m{d) and also 
hi and Pij (see Pen 2004) . In general it is not clear to what 
extend we can neglect e. What we do here instead is to 
consider cross-correlations. The galaxy-galaxy and galaxy- 
shear cross-correlations then relates directly to the matter 
auto correlations 



< Sg,5g^ > 

< ^aJs, > 

< Sg.Skj > 



hi < 5mi Srui > 

biPij < 5mi5mi > i < j 

bipij < Smi5mi > i < j 



(37) 



where we keep only leading order in Pij « 1. We have 
neglected < eiSi > and < ei€j > terms, which seems to be a 
good approximation on large linear scales as we flnd that the 
cross-correlation r defined in Eq |17| is very close to unity (see 
discussion in §2.4 and Figjl]). More generally, when r 7^ 1 we 
need to replace bi by hri in the second and thrid equations. 
This has little impact in our predictions as shown in j ]5.2.3[ 

Galaxy-galaxy cross-correlations from pairs of different 
redshift bins can be combined with the galaxy-galaxy auto- 
correlation to measure bias bi and pij and < Srm 5mi > with 
a radial precision given by the number of independent red- 
shift bins. The same information can also be obtained from 
combining galaxy-shear cross-correlation and galaxy-galaxy 
autocorrelation. 

Note that in our approach we ignore all the radial modes 
as the intrinsic < 5mi5mj > correlation is negligible for dis- 
joined top-hat bins i ^ j when Az > 0.02. It could be 
possible to include radial modes by using smaller redshift 
bins but this requires going beyond the Limber approxima- 
tion, a study that we leave for future analysis (see Challinor 
& Lewis 2011). In this sense our results are conservative, as 
we do not include radial modes or very fine radial bins. As 
we increase the radial resolution (or number of independent 
redshift bins) the number of independent transverse modes 
ki = £/ri also increases. 

Finally, recall that we do not include redshift space dis- 
tortions (RSD) in the modeling of angular correlations. This 
is a very good approximation for lensing because it has a 
very broad radial window which washes away the effect of ra- 
dial peculiar velocities. For the angular galaxy-galaxy auto- 
correlation in narrow redshift bins, this is not a good approx- 
imation (see Nock et al. 2010, Crocce, Cabre & Gaztanaga 
2010 and references therein). As mentioned in §2.3 our ap- 
proach is to model angular clustering as transverse modes in 
real space and include the effect of RSD through the ratios 
of the amplitude of clustering as we change from transverse 
to radial modes within a single redshift bin (see §2.5). 



3.2 Correlations &: Power spectrum 

Consider measurements of angular galaxy density and shear 
in a set of z-bins i = 1, ..,Nz. Then the projected measure- 
ment of A in bin i is: 



dz PAAz)6m{rO, z) 



(38) 
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where 5m is the 3D dark matter fluctuation, 9 gives the 
angular position, r — r{z) is the comoving transverse (or 
angular diameter) distance to redshift z and pAi is the ra- 
dial weight defined later. The label A refers to either the 
galaxy fluctuation, i.e. A — g, or a shear measurement of 
convergence, i.e. A — k. 

The observables to measure are the galaxy-galaxy (gg), 
galaxy-shear (gK) and shear-shear (kk) cross-correlations of 
angular fluctuations between redshift bins i and j, which are 
denoted by: 



(39) 
(40) 
(41) 



where 9 = \92 — 9i\ is the angular separation between 2 
measured fluctuations of galaxy density Sg or galaxy shear 
Sk (see Bernstein 2009 for generalizations of this). It is con- 
venient to express w{9) in terms of harmonic coefficients 

c{ey. 



w{9) 



C{£) 



2n 



d cos 9w{6) Li (cos 9) 



(42) 



(43) 



where Le{cos9) are Legendre polynomials of order £. Thus 
we have galaxy-galaxy, Cg^g-{1), shear-shear, CK.iK.j{£) and 
galaxy-shear, Cg^nj(l) cross power spectrum between bin 
i and j. Our galaxy sample is split into different galaxy 
types (or populations) that we label by gn, where n can 
take Nt values (we only use Nt = 2 here, but this can be 
easily extended), ie n = 1,2, ..Nt, for each population with 
different bias. We do not split shear by populations because 
we assume that shear is independent of bias. Thus the total 
number of observables for 2 populations, Nt = 2, is then: 

C.jii) = [C,l.,l^.(^),C,2,,2,(£),C,l.,2^.(£), 

C,i.«^.(^),Q2,.,(^),a,«,(^)] (44) 
For a generic cross-correlation we will use the notation: 

Ca^b.W (45) 

where A and B can take the values k for shear or gl, g2, ..gNt 
for different galaxy types. Note that while Ca^Aj = CajA^, 
in general: 

CA,s,(f) = Cg^A^e) (46) 
Ca^b^W / Ca.bM for Aj^B ; ij^j (47) 

In the small angle approximation (Limber 1954, Loverde 
& Afshordi 2008), there is a simple relation between C{£) 
and the 3D matter power spectrum P{k) =< 5'^{k) >: 



r{k,z) 



dz PA, (2) PBj (z) V{k, z) 



(48) 
(49) 



P{k,z) 
rH{z)r^{z) 

where rniz) = dr(z)/dz, and V is the a-dimensional power 
spectrum at k = [l -\- l/2)/r. In linear theory P[k,z) — 
D^{z)P(k). 



3.3 Lensing weights 

For shear or convergence fluctuations, i.e. A = k. the redshift 
distribution pA = Pk. in the above equation corresponds to 
the weak lensing efliciency between sources in redshift bin j 
and dark mater (lenses) at z: 



Pk^ {z) 



3QmHor{z) 
2H{z)a{z)ro 



dz 



I r{z';z) 
r{z') 



(50) 



where is the DM density at z = 0, ro = c/Hq, a{z) — 
1/(1 + z). Here r{z'; z) is the angular diameter distance be- 
tween z' and z and (j>j{z) is the normalized galaxy proba- 
bility density distribution for source galaxies in the redshift 
bin label j. For a uniform galaxy distribution in a top-hat 
window of total width A, : 



1/A, forl2-2j| < A,/2 



(51) 







otherwise. 



3.3.1 Narrow bin approximation 

For narrow redshift bins we approximate in Eq |50| by: 

P«.-P^,(-0-^f^'^^^^ ^or^<J (52) 
' 2H{Zi)ai rorj 

where = r{zi). This is an excellent approximation for the 
non-overlapping top-hat bins in our problem. For example 
for sources at Zi = 1, and a bin width of Azi = 0.1 (which is 
30 times larger than our default value at Zi = 1) the accuracy 
of this approximation is better than 0.7% for 2 < 0.95, which 
is all we need for non-overlapping z bins. This is illustrated in 
Fig|3) which compares the exact calculation (squares) with 
the approximation (lines across the squares) for Az^ = 0.05. 
This approximation greatly simplifles Eq |48| which becomes 
a single integral that can be done analytically for the galaxy- 
galaxy and the galaxy-shear case. Also note that p^ in this 
approximation is independent of the redshift bin width. 



3.4 Galaxy weight 

For galaxy type "n", i.e. A — gn, the redshift distribution 
Pa = Pgn in the Eq |48[ corresponds to the sum of the intrin- 
sic number density and the lensing magniflcation contribu- 



tion (see Eq 36 1: 

Pgrii (2) = b„. (2) -I- a„.pKi (2) 



(53) 



where 6„. = b„{zi) is the mean linear galaxy bias for galaxy 
type n at redshift bin i. The second term is due to weak 
lensing magnification with 

= a„{zi) = 5s„{zi) - 2 (54) 

where s„(2i) is the slope of the number counts Nn{< m; Zi) 
of galaxies of type n with apparent magnitude smaller than 
the survey flux limit m (plus whatever other cuts we do to 
the survey) at 2 = 2^ 

dlogioA/'„(< m;Zi) 



1(2,) 



dm 



(55) 
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Figure 3. Filled squares show the weak lensing efficiency pj.^ (z) 
in Eq |50| as a function of the lensed position z for a source at 
Zj = 1.0. The li ne a cross the squares corresponds to the ap- 
proximation in Eqj52| The other two lines show the weights p^^ 
and PKiPKj applied to the matter power spectrum V^k^z) in 
the shear-shear auto and cross-correlation of Eq |58| for Zi = 0.6. 
These distributions are always very broad, indicating that shear- 
shear can only measure the projected (2D) power. The shaded 
region shows the corresponding weights for galaxy-galaxy cross- 
correlation from MAG (or galaxy-shear), ie ipziPnj , for A; = 0.05. 
Here the power spectrum is measured in 3D with a resolution that 
is only limited by . 



Its value can be estimated using the same galaxy sample. We 
then have that modulates the weak lensing magnification 
effect of dark matter fluctuations at 2: < Zi. It is clear from 
this that both shear and magnification could measure the 
same information, modulo a„, noise and systematics, but in 
practice we need to study their correlations, which we do 
next. 



3.5 Galaxy-Galaxy (Magnification) 

For narrow and disjoint redshift bins we have from Eq |48| 
that for i <= j, 



+ 



'A, 



(56) 



where 5ij is the Kronecker delta and Vi = V{ki, Zi) in Eq 49 
with ki — {£ + l/2)/r{zi). When i — j the first term domi- 
nates (as Pa — 0), while for i 7^ j the second term is domi- 
nant (as 5ij = 0). The last term is always sub-dominant and 
will be neglected here|^ So we see here how galaxy cross- 



^ But note that if we can measure very well the two first terms, 
the last term contains all the shear-shear information, without 
need of direct shear measurements. 



correlations can be used to measure 3D power spectrum 
without the need of shear measurements. 

The nuisance biases b„;, can be measured by compar- 
ing dj with di. The main source of contamination here 
is intrinsic or induced correlations of galaxies in bin i with 
the ones in bin j. If these bins are well separated this can 
only occur because of photo-z transitions, ie galaxies moving 
from one redshift bin to another due to photo-z errors (see 



S3.ll and f3.12 for details 



The different pair combinations can be used to simul- 
taneously measure bi, Vi and Pij -To illustrate this claim, 
consider the ratio: dj /dk when i ^ j =^ k we have: 



dkW 



OlkPik 



(57) 



which are independent of bias and the power spectrum for 
any value of £. This is also true in the non-linear regime (i.e. 
with non-linear bias and non-linear dark matter clustering). 
This is therefore a direct geometrical measure. In a similar 
way one can build up ratios like Cu/dj to measure bi inde- 
pendent of Vi, or Cfj/di to measure Vi with independence 
of bias bi . 



3.6 Shear-Shear 

For the shear-shear term with i < j the cross-correlations 
are: 



dz PKiiz) PKj{z) V{k,z) 
dz l'3QmHo\^ r{zi;z)r{zj\z) 



(58) 
P{k,z) 



and C^iK ~ Cn-Ki for j < i. It is also possible to 
combine here different pairs of bins to get some tomographic 
(3D) information, but as illustrated in Figjs] the resulting 
weights are very broad, even if we use narrow galaxy bins. 
This means that the tomographic recovery is quite limited. 



3.7 Galaxy-Shear 

Finally, for narrow bins, we also find the galaxy-shear cross- 
correlation from Eq|48|for i < j: 



d 



bn^PijVi -\- a„^Cf, 



(59) 



and C, 



for i < j. Thus the off-diagonal information 
content of Cgmnj is identical (except for the global factor 
ctmi) to that of Cgmgmj, even when the noise and the sys- 
tematic could be quite different. 



3.8 3D P{k) recovery 

A comparison of Eq |59| or Eq |56| with Eq]58| illustrates how 
cross-correlations of galaxy-galaxy or galaxy-shear both re- 
cover the full 3D power spectrum, i.e. V{k), while the shear- 
shear only provides an integral over very broad redshift 
distribution. This is illustrated in Fig|3j which shows that 
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PKiPKj in Eq 58 is always very broad, while can be 



as narrow as needed. 

Another point to note is the fact that the galaxy-galaxy 
or galaxy-shear cross-correlations only depend onp^ ~ 10~^ 



while shear-shear depends on Pf. 



10" 



This means that 



the amplitude of shear-shear correlation is much smaller and 
is therefore sensitive to smaller amplitudes in systematic ef- 
fects. Moreover, the effect of systematics will tend to cancel 
in a cross-correlation analysis to a larger extend than in the 
auto-correlation. But these are generic considerations and 
ultimately the key point is that systematics are quite differ- 
ent in both measurements. 



3.9 Covariance 

In the Gaussian limit, the covariance between one pair of 
observables in Eq |44| at redshift bins (ij) and another pair 
at redshift bins (kl) is given by: 

e.,;,, = Cov(C.y,C,,) = (60) 

where A'^(^) is the number of modes at a given £. In our case 
we bin the I modes as A£ ~ 2//^^^ (see Cabre et al. 2007) 
to avoid correlation induced by the limited fraction of sky 
covered {fsky)- The number of modes of each £ is then 

N{£) = {2£+l)fskyA£~2{2£ + l) (61) 

The observables C in the covariance include observational 



total S/N is independent of the number of bins for large 
Nz- Something similar happens when considering RSD or 
in general 3D P{k) measurements. This means that there 
is no S/N gain in using many redshift bins. We could gain 
information in cases where the signal varies on radial scales 
comparable to the redshift bin width. Both cosmological pa- 
rameters and galaxy formation (i.e. biasing) vary on scales 
Aa — 0.1. This means that there is no information gain in 
using smaller bins in this case. 
In the case of galaxy-shear: 



^2 



N(£) N(£) 



(67) 



Here there are only galaxy-shear pairs per shear bin, 
so again the total signal to noise is quite insensitive to the 
number of bins. 

In the case of the autocorrelation C'u in Eq |65[ the total 
S/N increases with number of bins and there is a net S/N 
gain in using more bins. This gain is only limited by the 
increase of shot-noise for small redshift bins and the fact 
that very narrow nearby bins are no longer independent. 
Thus, while shear-shear and galaxy-shear are useful to re- 
cover galaxy bias, the gain from using narrow bins comes 
from galaxy-galaxy auto correlations. This is a key point to 
understand the results in this paper. 



3.11 Photo-z error transitions 



(62) 

(63) 
(64) 



where n„ and Tign are the surface density of galaxies with 
measured shear and galaxies of type n with good photometry 
respectively and is the variance in convergence from in- 
trinsic ellipticities. Note that (Jk < 1 while ng„ > so that 
the noise could be larger or smaller in shear-shear than in 
galaxy-galaxy, depending on the depth and quality of data. 



First we study the case of galaxy-galaxy cross-correlations. 
Consider the transition probability Tij that a galaxy at bin 
j is measured to be at bin i because of the photo-z error. 
The number of galaxies measured in bin i, Ni, is then: 



^3 Nj 



(68) 



where Nj is the true number of galaxies in bin j. To include 
the effect of photo-z errors in the cross-correlations we need 
to define the relative transition probability , or migration 
matrixrij , as the fraction of the galaxies assigned to z-bin i 
which really are in bin j, ie: 



3.10 Signal to Noise 

From the above covariance we can estimate the signal to 
noise S/N ratio for the galaxy auto-correlation Cg-g - {£) 



(5/iV)L. = 



N{£) 



(65) 



Til 



where by construction: 



X] = 1 



Vi 



(69) 



(70) 



and compare it to the galaxy-galaxy cross-correlation 

Cg,g,{£) 



(S/N) 



9i3j 



C?,N{£) 



oc N{£)AiAj oc 



N{e) 
m 



(66) 



where in the first step we have used that C^j < CjjCu and 
Cii oc 1/Ai, where cx l/N^ is the redshift bin width 
and Nz is the number of redshift bins. The total S/N is the 
sum over all the Nz redshift bins. In the case of the cross- 
correlation there are Nz{Nz — l)/2 pairs and therefore the 



Note that when rij is a smooth function of the redshift dis- 
tance Zi - 



rij will also be independent of the redshift 



width Az for Iz, 



> Az. These ra transitions cor- 



respond to the contamination matrix Cps in Bernstein & 
Huterer (2010) and give the probability P{zs\zp) for a true 
redshift is Zs to be measured in photo-z redshift Zp. The last 
equality in Eq |69| just indicates that the same probabilities 
apply for the mean density as for regions with fluctuations: 
Ni =< Ni > {1 + Si). We then have that the correlation 
Wij =< SiSj > in photo-z space can be related to the true 



© 2011 RAS, MNRAS 000, 1-10 



14 Enrique Gaztanaga, Martin Eriksen, Martin Crocce etal. 



correlation Wij =< SiSj > as (see also Benjamin et al. 2010) 

= <(l + ^0(l + ^i)>-l= t^'^'^ 1 

<Ni >< Nj > 

= ^^rikrji{l + Wki) — 1 ^^^Tikrji Wki (71) 

kl kl 

Thus, photo-z errors result in a mixing of the cross- 
correlation measurements: 



Cij = ^ rik rji Cki 



(72) 



where dj is the observed cross-correlation in photo-z space 
and Cki is the true cross-correlation. The new covariance 



matrix (i.e. Eq 60 1 will now be 



e = Cov{C\ C) 



(73) 



Note that noise (e.g. shot-noise) should be added to the ob- 
served (not the true) covariance. This could degrade the cos- 
mological parameters correspondingly. These considerations 
can be extended trivially to the case of the cross-correlation 
of 2 galaxy populations. We will then have different transi- 
tions rij for each population and assume no transition be- 
tween populations. In matrix notation: 



Cglg2 = rgl Cglg2 T 



(74) 



where rgi and rg2 represent the transition matrix for galaxy 
populations gl and g2 and T stands for the transpose ma- 
trix. In our case qi and 172 can take the values of Faint, gp, 
or Bright, gs, galaxies. 

For galaxy-shear and for shear-shear we use the same 
relations Eq |72|73l to account for photo-z contamination. In 
the case of shear, the observable is given by some mean 
over the number of galaxies where shapes are measured. So 
photo-z leakage has a different effect in shear fluctuations 
than in galaxy counts (see Bernstein & Huterer 2010). Nev- 
ertheless we note the shear correlations are also built from 
pairs of measurements and one can use the same relations 
as for galaxy-galaxy to account for photo-z contamination. 

Consider a photo-z distribution P(zp\z) giving the prob- 
ability that a galaxy at true redshift z is measured to be at 
Zp. For top-hat bins of width Ai and Aj we have: 



dz 

A,- J A, 



dzpP{zp\z) ~ / dzpP{zp\zj) (75) 

A,- 



where the integrals are around the corresponding top-hat 
bins centered in Zi and Zj. In the second equality above we 
have approximated the integral over bin j' by its mean value. 
Note that Tij 7^ Tji because the distribution P{z\zp) is in 
general different for different z (i.e. typically photo-z errors 
are a function of z). For a Gaussian photo-z distribution 
with error ~ (Jz{zj) 



2A, 

+ Ai 



A,/ 2 



erf{ 



2x + Ai 

'2z, 



erf 



erfi 



A. 



2a; - A, , 



(76) 



where Zij = zj 
(integration of a normal from to x) 



and erf{x) is the standard error function 



In general, it is possible to recover Vij from the ob- 
servables given some model for the intrinsic correlations, 
even without the combination of spectroscopic and photo- 
metric samples. The cross-correlation of spectroscopic and 
photometric samples provides another route to measure rij , 
extending the proposal by Newman (2008) to the cross- 
correlation in separate redshift bins. 

3.12 Uncertainties in the Photo-z transitions 

So far we have seem how it is possible to correct for photo-z 
errors. Here we will quantify how the uncertainties in our 
knowledge of the photo-z transitions Vij, which we will de- 
note by Aij, could affect our predictions. From Eq |72| we 
find that the error in C'a is: 



ACij — ^ ^ ( Raj Aia + Lie Aj 



(77) 



where we have defined the matrices R and L as: 

Lij = TikCkj ', Rij = "^^TjkCik (78) 

k k 

For a spectroscopic sample in the foreground (back- 
ground) we have that L ~ C (R ~ C). We will assume 
that this uncertainty is the same A,. — Aij for all pairs ij 
and that they are uncorrelated, so 

I 2 



< AijAki >= SikSjiAr 



(79) 



where Sik and Sjm are the Kronecker deltas. Note that in a 
more realistic situation errors could be different for different 
transitions, i.e. < Afj >^< A|j >, but here we just want to 
evaluate the overall size of this effect. One can easily extend 
this analysis to a more generic situation. The covariance 
caused by this uncertainty is then: 

Cov[af,Cki] =< Aa^ACki >= A^ (80) 



e 



ij;kl 



x[S,k(R^ R)ji + Sii{LR)kj +Sjk(LR)ii +Sji{LL^ 



evaluating the uncertainty in the Cij prediction, for a given 
value and error in Tij. This covariance is then added to the 
covariance matrix in Eq |73| that includes the observational 
noise, to explore how the FoM degrades as a function of the 
uncertainties in the bin transition A,.. This will give us an 
idea of how well we need to know rij for a given required 
FoM value. 

When we have different populations or cross-correlation 
of different observables (ie, galaxy with shear) we can easily 
generalize the above expressions, e.g.: 



ACa.s, = ( Ra^b, At + La^b^ a 



(81) 



where A'* is the uncertainty for type A. Types A and B 
can refer to density fiuctuations (e.g. A = gp or A = gB 
for Faint and Bright) and the same formalism can also be 
applied to shear. 



4 FIDUCIAL GALAXY SURVEYS 

We compare results for combinations of two galaxy samples: 
a deep {iAB < 24) sample, that we call Faint (F), and a 
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shallower (ms < 22.5) sample, call Bright (B). We consider 
four cases: the two samples alone and two ways to combine 
them. In one of the combinations (call F+B) we assume that 
the two samples are independent, ie over different parts of 
the sky. In the other combination (call FxB) we assume that 
they are over the same region of sky: typically one can think 
the B sample as a spectroscopic follow-up of the F sample. 

4.1 Faint Sample (F) 

The F sample is a photometric survey with a broad band 
photo-z error: 

22.5 < iAB < 24 and ~ 0.05(1 + z) 

similar to the upcoming DES (Dark Energy Survey, 
www.darkenergysurvey.org). We assume that we have very 
good photometry to be able to measure magnification 
and also shear (galaxy shape) information. But we will 
show results for magnification alone, shear alone and both 
combined. 

The cut at the bright end, iab > 22.5 is not an essential 
ingredient. It is included to make this sample separate in 
magnitude from the B sample, but all the results to follow 
are almost identical without such a cut. This selection results 



dN 
dfldz 



7.342 X 10*(^/0.467)i-^"e-(^/°-«") 



(82) 



over 0.1 ^ z < 2.0, where dN/dQdz is in units of number 
of galaxies per deg^ This is a fit to an evolving luminosity 
function that matches Blanton et al. (2003), Dahlen et al. 
(2005, 2007) and Jouvel et al. (2009). The above redshift 
distribution includes a completeness reduction of ~ 50%. 
For a survey of 200 deg^ it yields ~ 6 x 10^ galaxies. Note 
that we often refer to this sample briefiy as iAB < 24, but 
we always use a lower cut iAB > 22.5, so that galaxies are 
clearly separated from the Bright sample. This avoids poten- 
tial photo-z transitions between the two samples. The up- 
per line in Fig|4]illustrates the above model and compares it 
with the actual data in the public COSMOS photo-z sample 
(Ilbert et al.,2010) which has about 100000 galaxies over 2 
deg^ area with 22.5 < iAB < 24 and mean photo-z redshift 
of z = 0.9 for z > 0.2. Note that this area is quite small 
and subject to large sampling fluctuations (see Fosalba et 
al. 2008). 

For the fiducial bias in this sample we use four bias 
parameters 6f at Zi = 0.25,0.43,0.67 and 1.0 with linear 



interpolation to other redshifts (see S2.4I. The fiducial 
values for the 4 bias parameters follow: 

b[ = 1.2 + 0.4(2i -0.5), 

so that the bias approaches unity at z = and has a 
linear scaling, as shown in Fig[l] We then allow the 4 bias 
to vary without priors around these fiducial values. 

Note that we only consider WL and angular galaxy clus- 
tering probes for the Faint sample alone. It is also possible 
to measure RSD here even with poor photo-z (e.g. see Nock 
et al. 2010, Crocce, Cabre & Gaztanaga 2010, Ross et al 



8x101 



6x104 



-S 4x101 



Z 

T3 



2x101 




Figure 4. Upper and lower continuous lines show the number of 
galaxies per deg^ in our complete fiducial F (faint) sample with 
22.5 < lAB < 24 and the B (bright) sample with iAB < 22.5 
i.e. Eq |82| and Eq |83| multiplied by 2 to make then complete. His- 
tograms show the actual measurements in the COSMOS photo-z 
survey which has z = 0.90 and z = 0.57 respectively for z > 0.2. 
The long and the short dashed lines show the corresponding fits 
to iAB < 21.5 and iab < 23.0. 



2011, Crocce et al. 2011) but we have checked that the ad- 
ditional FoM is small compared to the one in spectroscopic 
surveys and decided to restrict to the more traditional use 
of photometric surveys. 



4.2 Bright Sample (B) 

The bright sample is a spectroscopic sample or very good 
photo-z without lensing information and defined by a flux 
limit: 

iAB < 22.5 
which results in: 



dN 
dQ,dz 



3.481 X 10*(2/0.702)^'0»^e-(^/°'^°^' 



(83) 



over 0.2 < z < 1.25 also including a completeness of 50%, 
as in the faint sample. For a survey of 200 deg^ it yields 
~ 2 X 10® galaxies. The lower solid line in Figlll illustrates 
the above model and compares it with the actual data in the 
public COSMOS photo-z sample (Ilbert et al.,2010) which 
has about 43000 galaxies with iAB < 22.5 with a mean 
photo-z redshift of z = 0.57 for z > 0.2. 

The fiducial photo-z errors will be CTz ~ 0.0035(1 -I- z), 
based on the PAU concept of using narrow band filters (see 
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§1). On linear scales, larger than 20 Mpc/h, this is almost 
equivalent to spectroscopic accuracy. In |5.8| we will explore 
the dependence of the results with photo-z error. Later on, 
in ij6j we will also consider a sample named "Bspec" with 
spectroscopic redshifts (a^ = 0) but with a lower density of 
1000 galaxy/deg'^. Dependence on density will be considered 
in section §5.6[ 

The four fiducial bias parameters, hf , for the B sample 
follow: 

&f = 2 + 2{z, - 0.5) 

Note that the bias at z; = is fe = 1 in both F and 
B samples. This is reasonable as bias only deviates signif- 
icantly from unity when the tracer is dominated by large 
masses or large luminosities, and this only occurs at high 
redshifts. 

For the bright sample we focus on RSD, which corre- 
sponds to the more traditional use of a spectroscopic sam- 
ple. BAO will also be considered in !|6] Using a good parent 
photometric catalog it would in principle be possible to do 
WL using shear or magnitudes over the same spectroscopic 
objects. The gain in such case comes from the combination 
of galaxy positions (with good radial information) and the 
weak lensing information in the parent sample. This belongs 
to the combination FxB below and is not considered as part 
of the features of B sample alone. The separation will allow 
for a more clear understanding of what we gain with the 
different combinations. 



case as this involves combining the redshifts information in 
the B sample with the WL information in the F sample over 
the same objects. This only happens when we combine both 
samples and the B sample is a subset of the F sample over 
the same area. This, together with the cross-correlation, are 
key and unique features of this combination. 

4.5 Redshift bins 



As discussed in §3.10[ the redshift bin widths have little 
impact in the RSD analysis or shear analysis, but they 
are important for the galaxy-galaxy auto correlations, i.e. 
CsiSi- We use non-overlapping top- hats with a width 4 
times larger than the photo-z error in the Bright sample, 
i.e. A = 0.014(1 + z). This corresponds to about 60 Mpc/h 
at z — 0.6 and is close to the smallest width where we can 
still neglect the intrinsic correlation between adjacent bins. 
For the Bright sample we will show that this width is large 
enough (as compared to photo-z error) so that we can al- 
most neglect photo-z errors. For the Faint sample, which 
has az = 0.05(1 + z), adjacent bins are strongly correlated 
because of photo-z errors. These correlations are taken into 
account in our analysis using photo-z error transitions (see 
S 3.11 1, so that our results are very similar to those obtained 



for larger bins, as expected. Thus the equivalent effective 
width is much larger for the F sample than for the B sam- 
ple, and is comparable to the photo-z error in the F sample. 
We will show in §5.7| how our conclusions depend on the 
redshift bin width. 



4.3 Independent Samples (F+B) 

Faint plus Bright: a direct combination (i.e. addition of 
Fisher Matrix) of the two cases above, assuming that both 
samples are independent and sample different parts of the 
sky. For a given area, the number of independent modes (or 
total information) available in this case is twice the num- 
ber of modes for the individual samples. In this respect one 
should consider that F-l-B has twice the area or cost as F or 
B. But given that we will be using different probes in F and 
B we do not quite have twice the information. 



4.4 Cross-correlation (FxB) 

Faint cross Bright: samples are combined over the same area. 
In this case the Bright sample is a subset of the original 
Faint sample, but here we consider disjoint magnitude bins 
to have independent sets of galaxies. This will avoid mixing 
of systematic effects, such as photo-z transitions. The total 
area sampled is somehow 1/2 of the case F-fB above, but we 
include the cross-correlations of the 2 samples for MAG (or 
WL) case. Here we use the MAG or WL and RSD probes 
in both the F and B samples as we have both adequate red- 
shift and photometric data . Traditionally one thinks of WL 
probes to be applied to deep photometric samples, such as 
sample F. As we will show below the combination of accurate 
redshifts and WL information in the B sample alone pro- 
duces figures of merit that can be comparable to the ones in 
the F sample. This information is only included in the FxB 



5 RESULTS 

In this section we will show the forecast on FoM produced 
by different probes on a particular choice of surveys. We 
will also explore how this forecast is affected by applying 
different assumptions on bias evolution models, photo-z er- 
rors, galaxy densities, survey areas and redshift bin-width. 
The goal is to understand the impact of each of the as- 
sumptions and to show how they affect the main results in 
this paper. We focus in MAG and RSD for most of this, 
but very similar results are found using other WL probes, 
as will also be shown. In i|6] we explore different surveys, 
including other WL probes (Shear-Shear, Galaxy-Shear) as 
well as BAO in addition to MAG and RSD, and compare 
expected constraints for different surveys. 



5.1 Bias and WL 

First we explore the relation between galaxy bias, the clus- 
tering from the galaxy-galaxy autocorrelation and the clus- 
tering from the different WL probes (i.e. shear and magni- 



fication). In [5.2 RSD is added to this picture. 



5.1.1 Bias fixed 

Fixing the galaxy bias is the same as assuming perfect 
knowledge of the bias. It represents the ideal case for re- 
covery of cosmological parameters given the galaxy samples. 
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bias: free-fixed 


Probe/Sample 


FoMji, FoM^ FoMu,-^ 


GG Faint (F) 
GG Bright (B) 
GS (B) 
GS+GG (B) 
WL-all (B) 


24-34 0.8-28 17-930 
27-53 0.1-31 4-1665 
8-11 0.04-6.2 0.3-70 
27-56 4-31 102-1749 
30-63 6-32 186-1998 



Table 2. Figures of Merit for DE equation of state wo — Wa 
(FoM„), growth 7 (FoM^) and joint (i.e. FoM„^). Each line 
shows different combinations of samples and probes: galaxy- 
galaxy auto-correlation (GG), galaxy-shear alone (GS) and 
their combination (excluding or including galaxy-galaxy cross- 
correlations). The first two entries correspond to GG over the 
Faint (F) and the Bright (B) samples. In each column we show 
two values separated by a dash "-" . The first value (bias free) uses 
a fit for b(z) using 4 free biasing parameters. The second assumes 
that bias is known (fixed bias). For 200 deg^ with Planck-|- SN-II 
priors. 



This case is not totally uninteresting as there are separate 
ways to constrain bias (i.e. from map reconstruction, halo 
modeling or other probes). It might also be possible to do 
direct lensing calibration of biasing (Bernstein & Cai 2011). 

Results for 200 deg^ with Planck and SN-II priors are 
shown as the right value in the pair of numbers separated by 
a dash (-) in Table |2] The first two rows compare galaxy- 
galaxy autocorrelations (GG) for the Faint (F) and Bright 
(B) samples. The Bright sample gives about a factor of 2 
larger combined FoM than the Faint sample, despite having 
2.5 less volume and fewer galaxies. This is because of the 
higher radial resolution and illustrates a key point in our 
approach: better radial accuracy of the foreground galax- 
ies improve both growth and cosmic history reconstruction. 



Bellow we show {\ 5.7 1 how this changes with the number of 
bins (or the bin width). 

The third entry is the results for galaxy-shear (GS) in 
the Bright sample. Here the combined FoM„^ is a factor of 
24 times lower than using the galaxy-galaxy autocorrelation 
over the same (B) sample. This shows the problem with 
weak-lensing: that has a broad kernel and is intrinsically 2D. 
Having good radial resolution in the foreground galaxies of 
the galaxy-shear cross-correlations does not help as much as 
in the galaxy-galaxy case (e.g. see ^3.10[ ). 

The next entry shows the combination of galaxy-shear 
with galaxy-galaxy autocorrelation, which is still dominated 
by the galajcy-galaxy results, but note how there is an in- 
crease in the FoM^ due to the geometric ratios pij which 
are only measured in galaxy-shear. The last entry include all 
weak lensing and galaxy cross-correlations, including mag- 
nification (we will call this WL-all). Again here we find a 
significant improvement, but results seem dominated by the 
galaxy-galaxy auto-correlation. Note how there is very lit- 
tle improvement in the FoM-y in the last two entries with 
respect to the galaxy-galaxy auto-correlation. When bias is 
known, weak lensing does not add new information about 
the growth over what we can measure with galaxy cluster- 
ing alone. 





bias: free-fixed 


Probe/Sample 


FoMji, 


FoM^ 


FoM„^ 


MAG (F) 


25-40 


1.3-30 


32-1198 


RSD (B) 


5-29 


6-16 


30-463 


MAG+RSD F-(-B 


28-138 


7-35 


180-4859 




MAG-f RSD FxB 


86-197 


21-39 


1760-7733 





Table 3. Same as Tablc[2] for magnification (MAG) and Red- 
shift Space Distortions (RSD). The first two entries correspond 
to MAG over the Faint (F) and RSD over the Bright (B) sam- 
ples. In the third row (F-l-B) the MAG and RSD probes over the 
F and B samples are combined assuming to be uncorrelated and 
to come from different regions of the sky. In the last line (FxB) 
the two samples are combined as being over the same area of the 
sky, including its cross-correlation and covariance. 



5.1.2 Bias free 

This picture changes dramatically when we allow bias to 
vary. The bold value on the left (in the pair separated by a 
dash) in Table |2] shows how the above results change when 
we measure the bias evolution 6(z) without priors. As justi- 
fied in §2.4| we use 4 free bias parameters for each population 
(F or B) spaced by Aa ~ 0.1 (i.e. At ~ IGj/r) and interpo- 
late linearly in redshift between these 4 points. 

In general, there is significant reduction in all FoMs. 
Notice how FoM^ is more affected than FoMi„ because bias 
is degenerate with growth in the power spectrum amplitude. 
In the case of Galaxy-Shear the FoMi^, is hardly affected by 
the bias. Recall from Eq |59| that in this case the angular 
power spectrum between galaxies at redshift bin i and shear 
at J > i is Cg^Kj ~ biPijVi. The information about 'w{z) is 
encoded in pij but is also degenerate with the growth infor- 
mation which depends on w{z) through Q{z). This produces 
a degeneracy in 'w{z) for GS that is not broken by knowning 
the bias. These effects are less severe for GG because of the 
better radial resolution which helps breaking this degener- 
acy even when bias is free. Note how in this case FoM^,^ is 
now better for the Faint sample than for the Bright sample 
in contrast to the situation when bias is known. 

Note here how the WL-all combination (last entry) pro- 
vides a significant improvement in FoM^ with respect to the 
GG (B) or GS-I-GG cases. In contrast, when bias is known, 
there is no improvement. This improvement comes from 
bias, which can be measured with the combination of galaxy- 
shear and galaxy-galaxy statistics, eg bi oc Cg-g^/Cg^K . This 
is a key point to notice. So far these results do not include 
the cross-correlations between the Faint and Bright samples 
which better exploit the complementarity of deep weak lens- 
ing measurements with good radial resolution in the fore- 
ground lenses. They also do not include RSD which provides 
another way to measure bias. We will explore this in the fol- 
lowing section about magnification (MAG) and RSD and 
also later on (^ when we combine all the different probes. 



5.2 Combining RSD and MAG 

Table [3] shows a comparison of the different FoM for MAG 
and RSD probes. As above, we will focus first on the results 
for a fixed bias. 
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^0 7 

Figure 5. Contours (Ax^ = 2.3) of too ~ — 7 in magnification (MAG) over a Faint sample with 22.5 < iab < 24 (labeled F200, 
blue) as compared with redshift space distortions (RSD) over a bright subset (labeled B200, purple) with < 22.5 with very good 
photo-z resolution: = 0.0035(1 + z). Bias evolution b{z) and other 6 cosmological parameters are marginalized over. The red contours 
correspond to the combination of both probes over different sky areas (F+B), while the yellow contours show the combination over the 
same area (FxB), including also all cross-correlations between galaxies in both samples. 



5.2.1 Bias fixed 

As shown in Table [Sj MAG with fixed bias produces better 
results than RSD. The combined FoMu,.y is a factor of 2.6 
times better for the Faint sample using MAG. 

One of the main points we want to stress in this paper 
is shown when we combine MAG+RSD. We can see a large 
improvement both for the separate samples (F+B) and for 
the cross-correlation analysis (FxB). This is highlighted by 
a box in Table |3] Note how the improvement is important 
for all FoMs. 

The F-l-B combination improves the FoM^^ and FoM„ 
by factors of 3-8 and 3-4, while the gain in F0M7 is more 
modest. These results will change dramatically when we al- 
low the bias to be free. 

For the F-l-B combination with fixed bias, having two 
separate samples over different regions of the sky is worse 
than doing the analysis over the same region, even when 
we are effectively doubling the area in the first case. The 
difference between F-l-B and FxB is almost a factor of 2 for 
fixed bias. This is because FxB includes more observables 
than F-fB. As explained in Section 4.4, in the FxB case we 
use the MAG and RSD probes in both the F and B samples 
and we also include the cross-correlation of the F and B 
samples. For F-l-B we only consider RSD in B, MAG on F 
and no cross-correlations. 



5.2.2 Bias free 

The above picture changes when we allow the bias to be free. 
In all cases the FoM degrades significantly. The degradation 
depends on the sample, probe and FoM considered, as il- 



lustrated by bold numbers in Table |3] In MAG the FoM^ 
degrades more than FoM^,, while the opposite happens in 
RSD. Note how fixing bias improves FoM„ for RSD (B). 
This is because of the improvement in D{z) which allows 
w{z)to be measure through its dependence on Q,{z). This is 
in contrast with GS in Table [2] which shows little improve- 
ment of FoMu, when bias is fixed because of further degen- 
eracies with the lensing geometrical parameters in pij. The 
combination of MAG and RSD can break these additional 
degeneracies. 

When we correlate samples over the same region of the 
sky (FxB) we find much better results than when the sam- 
ples are independent (F-l-B). The improvement is about a 
factor of 3 for FoM™ or FoM-y and 10 for FoM^-y. Bias is 
better constrained when both samples are over the same re- 
gion as we can include the cross-correlations between then 
and use both MAG and RSD probes in both samples. So 
even when there is less information in terms of number of 
independent modes or total area in FxB than in F-l-B, there 
is more constraining power in the FoM of cosmological pa- 
rameters. This is illustrated in Fig[5] which shows the 3 pro- 
jections (wojWa.y) in the FoM^,^ ellipsoid. Note how in the 
Wo — Wa plane there is little gain in the MAG+RSD (F+B) 
combination as the improvement is in 7. In FxB the de- 
generacies are strongly reduced thanks to the bias measure- 
ments. 

For the F+B comparison there is also a significant im- 
provement with respect to the separate F and B cases. In 
this case the relative improvement is smaller when bias is 
free (as compared to the case of known bias). It comes from 
having different biases for the F and B samples, but no cross 
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measurement to relate them, as in the FxB case where we 
use MAG cross-correlation to relate the bias in F and B. 

Allowing the bias to be free reduces the combined 
FoM„^ by factors of 37 and 15 in MAG and RSD. But this 
reduction is only a factor of 4.4 for the FxB combination. 
These results assume no priors on bias. We will check in next 
sub-section how this changes when we have some priors on 
bias. 



5.2.3 Bias stochasticity 

In section §2.4 and Figjl] we have argued for r being very 
close to unity in numerical simulation. Therefore through 
the paper, the bias is assumed deterministic (r=l), but here 
we in addition test how a stochastic bias would impact our 
result. A stochastic bias is introduced by replacing Eq|56|by 



(84) 



where r„. is the galaxy-matter cross-correlation coecient for 
the Hi galaxy population. The variation of r„. with redhift 
is parametriced in the same manner as the bias 6„;, using 
4 parameters which are scale independent on linear scales. 
We use = 1 as the fiducial values and allow r„; to vary 
without priors. 

When marginalizing over these new parameters, FoM„^ 
changes by less than 4% for the Faint sample.. The impact 
on the Bright sample is even smaller and under 1%. This 
small effect can be understood by comparison to the GG 
result for the F sample in top right corner of Table |2] most 
of the FoM in MAG comes from the GG contribution (the 
first of the above equations), which is independent of r„.. 
Such a small dependence on r„. will have little impact on 
the measurement of bias or on the FxB cases. 



5.3 Priors on bias 

Figure [6] shows how the FoMjn-/ increases as we reduce the 
priors in the bias parameters. These priors can come from 



other observations or from modeling of bias (e.g. see S2.4|. 
The same gaussian prior is used for all the bias parameters. 
F-l-B adds MAG in the F sample with RSD in the B sample, 
while FxB include all combinations and cross-correlations. 
The figure shows that the FxB combination is less affected 
by priors than the other combinations and reflects the fact 
that bias is better measured by the FxB combination. A 
prior of 0.01 in bias improves the FoM of FxB and F+B by 
about a factor of 2.4 and 3.8 respectively for the 5000 deg^ 
sample. With such priors, the FxB combination still has 
FoMu,7 that is 3.2 times larger than F-l-B. This is because 
in for the FxB combination we use both RSD and WL over 
the bright and faint sample, while in F-|-B we only use RSD 
over B and WL over F. 

Also note how for no priors the FoM in the small sam- 
ple combination FxB200 is almost as good as the F-I-B5000 
combination which has 25 times larger area. 

If we use 2 biasing parameters per population instead of 
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0.01 0.1 

Prior in bias 



Figure 6. FoM,i,^ in the MAG-I-RSD combination as a function 
of priors in the bias parameters (absolute rms error). Dashed lines 
correspond to the FoM when we combine samples from different 
regions of the sky (F-l-B), while continuous lines correspond to 
the cross-correlation of samples over the same area (FxB). Up- 
per/lower lines corresponds to 5000/200 deg^ surveys. 
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Figure 7. Variation of the different FoM: FoM^ (top line), 

FoM'^^j and FoM^, (bottom line) as a function of the Survey 
Area. This is RSD-I-MAG combined over the same area (FxB). 
The dashed line shows the scaling with area for reference. 



4 to describe b{z) (i.e. see [ 2.4 1 we find an increase of about 
40% in the FxB combined FoM„^ for both 200 or 5000 deg^ 
samples. 

We conclude that the impact of bias evolution is well 
understood and our results are quite robust with respect 
bias modeling. 
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Figure 8. Lines sliow tiie combined FoM,„^ as a function of tlie 
magnitude limit of the bright sample for different samples and 
probes, as labeled. The top line corresponds to the FxB cross- 
correlation using both MAG and RSD. The lower lines correspond 
to MAG alone over the Bright sample with 5000 deg^ or 200 deg^ . 



Figure 9. FoM^-y as a function of the relative galaxy density n 
in B200 with respect to the fiducial values (in fj4| for different 
probes, as labeled. As usual, F+B adds MAG in the F200 sample 
with RSD in the B200 sample, while FxB include all combinations 
and cross-correlations. The dashed lines show a fit to fi^. 



5.4 Survey Area 

Figure [7] shows how the MAG-I-RSD combination over the 
same area (FxB) changes as a function of the survey area. 
For the Fisher Matrix without priors we expect the follow- 
ing scahng: FoM^, oc A, FoM^. oc A'^^^ and FoM^,^ oc A^^^, 
according to Eq |14[ as it depends on the number of param- 
eters in each FoM definition. The values shown in Fig[7]are 
rescaled accordingly and show the FoM corresponding to 
the geometrical mean of the parameters involved. After this 
rescaling, they should all scale linearly with the area. But we 
are using Planck and SN-II priors which tend to dominate 
the results for smaller area, resulting in a weaker dependence 
in A. A fit to Fig[7] yields: 

FoM^ cx A°-^'^ (85) 
FoM., cx (86) 
FoM^., cx A°-^^ (87) 

Note that FoM^ is larger than FoM^,, which indicates 
that 7 is better determined than the geometrical mean of 
Wo and Wa- This is despite the fact that we include Planck 
and SN-II priors, which include priors in wo and Wa, but no 
priors in 7. 



5.5 Bright Magnitude limit 

Figure [S] shows how the FoM changes as we change the mag- 
nitude limit (mi) of the bright sample. For A''(2:) we use the 
distributions shown in Fig|4] The combined FoM^^ for MAG 
in the Bright sample depends strongly on mi, as shown by 
the bottom lines in the figure, for both 200 and 5000 deg^. 
When combined with the Faint sample and RSD, the depen- 



dence is weaker. For the latter case we find: 

FoM^^ cx L26'"' (88) 



5.6 Shot-Noise 

What happens if we reduce or increase the galaxy density? 
The bright sample {mi < 22.5), for which we so far have 
assumed a PAU-like survey measuring low resolution spectra 
for most of the galaxies, likely has a lower density for a truly 
spectroscopic survey. A different density can also come from 
only considering a subset of the data or dividing the galaxies 
in types. 

Figure |9] shows how the FoM from the different probes 
changes as we reduce (or increase) the density in the bright 
sample with 200 deg^ (B200). The dilution is given in terms 
of the default density in Q which already corresponded 
to a 50% completeness. Both RSD and MAG are quite af- 
fected by sample dilution, especially when we approach val- 
ues lower than 10% completeness. The data can be fitted 
with a power-law in the relative density n'', with (3 ~ 0.15 
for MAG and /3 ~ 0.3 for RSD (dashed lines in the figure). 
It is therefore important to keep a high density if we want 
a high FoM, specially for RSD. 

The effect is less severe in the F-I-B200 combination 
{P ~ 0.15) because the density of the faint sample (F200) 
is kept constant. For the FxB200 combination the effect is 
in between the one for B200 and the one for F-I-B200. The 
effect of dilution in the FxB200 case can be well described 
as: 

FoM^'f cx n°-2 (89) 
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Figure 10. Combined FoM™^ as a function of the relative red- 
shift bin width (with respect to the fiducial value of 0.014(l-|-z)) 
for the different samples and probes. As usual, F+B adds MAG 
in the F sample with RSD in the B sample, while FxB include all 
combinations and cross-correlations. 



5.7 Redshift bin width 

Figure [To] shows how the FoM^n, changes when we change 
the redshift bin width from 4 to 20 times the photo-z error 
in the Bright sample, i.e. 1 to 5 times the fiducial values. 

Changing the bin width has little effect on the Faint 
sample for MAG. The fiducial z-bin, 0.014(1 + z), is al- 
ready smaller than the photo-z error in the Faint sample, 
0.05(1 -I- z), indicating that nearby bins are strongly corre- 
lated. Increasing the bin width by a factor of 4 only changes 
FoMu,^ by 12%. This illustrates that the method to account 
for transition probabilities in §3.11| and their covariance in 
§3.9[ properly take into account the impact of photo-z errors 
in the MAG forecast. We do expect some small degradation 
as we increase the redshift bin width because of the lost of 
radial information. This degradation is more significant for 
a bin width which is larger than the photo-z error. 

For the Bright sample different redshift bins are inde- 
pendent as they are at least 4 times larger than photo-z 
error. The impact of changing the bin widths is quite differ- 
ent for RSD and MAG. For RSD the results are independent 
of the bin size, once there are enough bins. This is because 
we are using 3D information within each bin to measure 
the distortions: this is possible because we are in the linear 
regime and we assume photo-z accuracy that is good enough 
for this task, see §2.5| The redshift bin slicing only affects 
the radial resolution in the parameter estimation, but it does 
not add new modes as happens when we do 2D clustering. 
So as long as there are enough bins to capture the slow cos- 
mic evolution the results are unchanged in 3D. The F+B 
case suffers little degradation because it is just the sum of F 
with MAG and B with RSD, neither of which changes with 
number of bins. 
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Figure 11. Symbols show the combined FoM^-y as a function 
of the relative photo-z of the bright sample (with respect to the 
fiducial values in <Tz = 0.0035(1 -I- z)) for the different samples 
and probes, as labeled. Dashed lines show the best fit to e~^"^ . 



For MAG in the Bright sample the effect is more im- 
portant because the result is dominated by the 2D galaxy- 
galaxy (GG) correlations (see Table[2|, whose signal-to-noise 
depends linearly in the number of redshift bins (see i 3.10 1. 
If we reduce the number of independent bins, we directly 



reduce the FoM. In the case of MAG for the Bright sample 
a bin width that is twice as large reduces the FoMm^, by a 
factor of two. This agrees with the reduction of the 2D GG 
FoM for the B sample shown in Fig. |10| The degradation is 
smaller in the GS because its signal-to-noise is independent 



of the number of bins (see S 3.10 1 



The influence of the 2D galaxy-galaxy contribution is 
also apparent in the top lines of Fig. |10[ which shows the 
FxB combination for MAG or for WL-all (which includes all 
galaxy and Shear correlations). As we decrease the redshift 
bin width by a factor of 4, the combine FoM of FxB in 
MAG+RSD increases by 75%. 

Overall, this illustrates well one of the main points we 
want to stress in this paper: that the benefit of the FxB com- 
bination comes from the combination of many independent 
redshift bins. 



5.8 Photo-z errors 

Figure [iT] shows how the combined FoM^,^ changes as we 
change the relative photo-z errors, a^, in the bright sample. 
Here the photo-z error distribution is Gaussian. 

The variation on FoMm^ is smaller for MAG than for 
RSD. This makes sense as the bright MAG analysis uses rel- 
atively broad redshift bins: Az = 0.014(1 + z), i.e. 4 times 
the nominal photo-z error of (Tz = 0.0035(1 + z). This makes 
our result conservative but quite insensitive to photo-z er- 
rors. In the case of RSD the degradation comes from the 
S/N reduction on modes affected by the photo-z error. 

For the FxB results the effect is less severe as here we 
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Figure 12. Combined FoM^^ as a function of the absolute error in our knowledge of the photo-z transitions for 5000 deg^ survey. 
Left: Results correspond to cosmic magnification (MAG with free bias, continuous lines, or fixed bias, dotted line) and Shear-Shear 
cross-correlation (Shear, dashed lines) alone, i.e. without RSD, for both the Faint (F) and Bright (B) samples. Right: Results correspond 
to the FxB cross-correlation of the F and B sample, including RSD with different WL probes: magnification (MAG with free and fixed 
bias, top and bottom continuous line), Galaxy-Shear cross-correlations (GS, short-dashed), Shear-Shear cross-correlation (Shear, dotted 
line) and all WL combined (long dashed lines). 



also include results for the faint sample for which the photo-z 
error is fixed. These tendencies can be fitted to an exponen- 
tial decrease: e~'^°'^ . For both the 200 and 5000 deg'^ we find 
that to a good approximation: 



(90) 
(91) 



FoM^^^« oc e-o.ik./o.oo35(i+.)]^ 

These fits are shown as dashed lines in the figure. MAG in 
the B200 sample falls a bit less steeply {A ~ 0.07) than 
MAG in the B5000 {A ~ 0.1). 

Note in Figure [TT] how, despite having 25 less area, the 
combined FxB200 result is better that the separate MAG 
or RSD results over B5000. This is also true for results over 
the F5000 sample. 

5.8.1 Uncertainties in photo-z transitions 

So far we have considered the case where we know perfectly 
well how photo-z uncertainties change the observables. This 
corresponds to the case Ar = 0, where Ar is our uncertainty 
in the fraction of galaxies rtj that are wrongly allocated to 
a given redshift bin, as defined by Eq]79| The left panel in 



Fig 12 shows how the FoM for Magnification (MAG) and 



Shear-Shear cross-correlations (Shear) degrades as we in- 
crease the absolute error Ar. In the left panel we only in- 
clude WL probes and no RSD. 

For MAG, the degradation in the FoM is similar for the 
Bright and Faint samples. In the limit of small errors in r^j 
the FoM of MAG is a 5.6 times larger for the Faint than for 
the Bright sample (bias free), because the former is much 



deeper, has more redshift range, volume and galaxies. Thus, 
the better radial resolution in the Bright sample does not 
quite compensate for its smaller volume when we also have 
to fit for bias. But when bias is known (dotted lines) the 
Bright sample produce 2 times better results, which indi- 
cates that the larger volume in F is fully compensated with 
better radial resolution in B.. As we increase the uncertain- 
ties in the photo-z transitions the FoM degrades quickly, as 
we loose the capacity to locate radial bins. 

Recall that the Faint sample has a much larger photo-z 
error than the Bright sample and therefore the fraction of 
galaxies rij will also be much larger. This means that for 
a fixed absolute error in rij the relative error will be much 
smaller for the Faint sample. Does it make sense then that 
the FoM is degraded equally for the Fain and Bright sample, 
as shown in Fig |12| ? 

As discussed in §5.7| results for MAG in the Faint sample 
depends very weekly on the redshift bin width as long as it 
is not much larger than the photo-z error. To understand 
our new results we can increase the bin width in the Faint 
sample to make it similar or slightly larger than its photo- 
z error. In this case we have that for both the Faint and 
Bright sample ru ~ 1 and rij < 1 for i ^ j. In this case, the 
variance in dj caused by an absolute error Ar in rij is: 



{Aa,fc^AUcl + C^j) 



(92) 



This should be compared to the variance in dj in Eq |60| In 
the limit Ca ~ Cjj > dj we have: 



{Adjf 

Cov[dr,c^^ 



2AIn{1) 



(93) 
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which is independent of dj and therefore of az[^If we want 
this ratio to be smaller than 1/4, so that the error in dj 
is not dominated by Ar, we find that we need Ar < 10~^ 
for £ ~ 300, which is in good agreement with Fig |12| and 
illustrates why this result is the same for the F and B sam- 
ple, and why its shape is also independent of the redshift 
bin width. We have checked this prediction by redoing the 
results on the left panel of Fig |12| for twice the bin width. 
For the Faint sample we find a very similar result with a 
FoMu,^ that is only 5% lower for the larger bin width for all 
values of A^. For the Bright sample the results are a factor 
of 2 lower, as expected from §5.7[ with a similar shape as 
function of Ar . 

For absolute errors larger than about 1%, the degra- 
dation for the shear-shear cross-correlation (dashed line la- 



beled Shear in left panel of Fig 12 1 is smaller than for MAG 



This is because MAG provides better constraint than Shear 
alone (when bias is fixed). When bias is known, MAG has 
a larger FoM than Shear for Ar < 1% because of it's bet- 
ter radial resolution. This is specially true for the Bright 
sample, which also has too small a volume to take advan- 
tage of shear. But because of the broad radial resolution, 
the FoM from Shear is less sensitive to the uncertainties in 
the photo-z transitions. Note that this formalism for degra- 
dation is more general (for narrow bins) than the one based 
on the uncertainties in the centroid position of the photo-z 
source distribution (eg Huterer et al 2006). 

For the MAG probe, we need to know the transitions 
to better than 4% (1%) if we want the FoM to degrade 
less than 50% (11%). This is the case for both the Faint and 
Bright 200 deg^ samples. For the 5000 deg'^ samples we need 
to know the transitions to better than about 2.5% (1%) if 
we want the FoM of MAG-I-RSD to degrade less than 50% 
(13%). For WL-I-RSD a 1% uncertainty results in a 50% 
degradation. A bright sample with 5000deg^ where we only 
know transitions to 50% gives similar FoM in MAG as a 
sample with only 200deg'^ where transitions are known to 
be better than 2%. 

The relative impact in the FoM is similar if we con- 
sider combined probes with RSD, which we have assumed 
is not affected by uncertainties in the photo-z transitions. 
This is illustrated in the right panel of Fig |12[ which shows 
the corresponding results for the FxB combination of RSD 
with different WL probes. Here we take Ar to be the same 
for the F and B sample, which is probably pessimistic for 
bright spectroscopic samples (but not for PAU-like sur- 
veys). The short-dashed line corresponds to the Galaxy- 
Shear (GS-I-GG) cross-correlation (which also includes the 
galaxy-galaxy autocorrelation Cg-g-) while the continuous 
lines correspond to MAG (i.e. Cg-g ). The dotted line shows 
the shear-shear results. The long dashed line combines (with 
appropriate covariance) all WL lensing probes (ie MAG, GS, 
GG and SS). As can be seen in the figure all probes are af- 
fected in similar way by the uncertainty in the photo-z tran- 



^ If we include shot-noise we have to divide this by a factor [1 -|- 
l/{nCii)]^ which is typically of order unity, but depends weakly 



Probe/Sample 


FoM„ FoM^ FoMu,^ 


FxB no wiggles 
FxB wiggles 

F-f B no wiggles 
F-|-B wiggles 


78 (5.6) 21 (6.3) 1595 (16) 
86 (2.8) 20 (5.5) 1760 (35) 
26 (0.4) 6.5 (3.2) 170 (0.6) 
28 (0.2) 6.5 (3.2) 180 (1.2) 



Table 4. Impact of having BAO wiggles in the MAG-I-RSD com- 
bination for 200deg-^ survey with Planck and SN-II priors. Results 
without priors are in parentheses. 



sitions. When the uncertainty in the bin transition Ar are 
less than 0.1 we find: 



FoM^j^ oc e 



(Ar/0.18){A/200)" 



(94) 

where A is the survey area in deg^. The A dependence re- 
flects the fact that to achieve a given FoM degradation we 
need to know the transitions around 15% better (i.e. Ar 
should be 15% lower) for A=5000 deg^ than for A=200 deg^. 
This makes sense, as the errors are lower for A=5000 deg'^, so 
the requirements should be higher. Accuracies of few percent 
in Tij seem within reach with current or near future photo-z 
codes with appropriate calibration, so this is not a critical 
limitation to measure magnification. 

Note how the FoM can increases by over a factor of 100 
when we cross-correlate the samples and include RSD, i.e. 
compare the values of the FoM of FxB in the right panel with 
the FoM of F or B separately in the left panel. Also note 
how magnification becomes better than shear when com- 
bined with RSD (ie compare left and right panels). 



5.8.2 Non-Gaussian transitions 

We have also explored the case of non-Gaussian transitions 
by running a photo-z code over galaxy simulations and ex- 
ploring how this change the results with respect to the Gaus- 
sian case. We have estimated the photo-z transitions rij us- 
ing the photo-z code over mock galaxy simulations. Details 
of this study will be presented elsewhere. Here we just want 
to point out that in many cases we find that surprisingly the 
FoM increases for the Non-Gaussian transitions compared to 
the Gaussian case. The reason for this is that for the Gaus- 
sian comparison we use a gaussian rms az corresponding 
to the the 68% confident level of the photo-z distribution 
in simulations. The actual (non-Gaussian) photo-z distribu- 
tion is often much sharper in the center and has longer tails 
than a Gaussian. If the distribution is known, this results 
in a better radial resolution and therefore better FoM. In 
reality this improvement is limited by our uncertainties in 



the photo-z distribution (see S3.12 and below) 



on the sample and redshift bin width. 



5.9 Impact of BAO wiggles 

The BAO wiggles in P{k) can be measured in the angu- 
lar galaxy-galaxy correlations when we use narrow redshift 
bins. We test this by comparing results using the standard 
Eisenstein & Hu (1998) P{k) with and without the BAO 
feature for the same cosmology. 

Results for MAG+RSD are shown in Table g] There is 
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Parameter 


MAG 


RSD 


MAG+RSD 




21.1 


76.2 


27.1 




48.7 


69.2 


50.9 




19.3 


75.8 


18.6 


h 


19.8 


76.0 


21.2 


ris 


0.4 


0.0 


40.1 


o-g 


0.0 


0.9 


1.1 


&f 


77.6 


28.4 


2.6 


b? 


79.1 


31.9 


3.0 




81.5 


40.5 


5.0 


»§ 


82.0 


46.0 


5.6 


bi 


92.0 


66.7 


47.3 


bi 


92.0 


67.7 


46.1 


bf 


93.9 


78.3 


62.8 


bf 


94.0 


79.7 


64.5 



Table 5. Degradation (per cent reduction) of the FoM^-y when 
we add a new parameter to the FxB200 cross-correlation forecast. 
The change is relative to the case where we have all parameters 
varying except for the one quoted (without priors). The 4 bias 
parameters for the Faint and Bright populations bp and bg are 
the parameters that gain most in the MAG+RSD combination. 

no impact on RSD as the shape is fixed in tliis case. Without 
priors (in parentheses), the combined FoM in MAG+RSD 
(FxB) increases by a factor of 2 due to this effect. Of course 
this change mostly affects w{z) FoM as 7 is measured by the 
amplitude and not the shape of P{k). We find very similar 
relative improvement in FoMu, for SOOOdeg^ without priors. 
Note how the F+B combination also increases by a factor 
of 2 with BAO wiggles. This is despite that in F+B we 
use MAG for the faint sample alone, which has low radial 
resolution, which translates into a lower resolution in k for 
the shape measurement. 

The increase in FoM„ with priors is only 10% for FxB 
for 200deg^ and 31% for SOOOdeg^. Even with priors, these 
are quite substantial gains and illustrate that the BAO wig- 
gles are indeed well measured with these probes. 

5.10 Degeneracies 

We investigate here the gain in the MAG+RSD (FxB) com- 
bination for different parameters. Table |5] shows how the 
FoMw-y (with Planck and SN-II priors) degrades when we 
add a new parameter. We see how the bias parameters are 
the ones where the MAG+RSD combination degrades less 
than the separate RSD and MAG. So they are the ones that 
benefit most from the MAG+RSD combination. For exam- 
ple note how allowing 61, the bias at the lower redshift of 
Zi = 0.25, to vary causes the FoM to degrade by 30-80% 
for RSD or MAG separately but only by 3% for the FxB 
combination. 

In some cases, like Ua, the degradation increases when 
we consider the combined RSD+MAG probes. This is be- 
cause the result for separate probes is dominated by the 
priors. In this case adding a new parameter does not de- 
grade the FoM. The combination of probes provide better 
constraints which are no longer dominated by priors. This 
results in an increase in degradation. In the case of ag, the 
degradation is always low as its value is dominated by priors 
and does not affect much our FoM. 





MAG 


RSD 


MAG+RSD 




F200 


B200 


F+B200 


FxB200 


FxBSOOO 




5.3% 


47% 


5.5% 


3.3% 


1.5% 




1.8% 


9.2% 


1.8% 


1.2% 


0.6% 


h 


2.8% 


24% 


2.8% 


1.8% 


0.83% 


OS 


0.9% 


0.9% 


0.9% 


0.9% 


0.76% 




5.8% 


47% 


5.7% 


3.6% 


1.88% 


wo 


15% 


17% 


14% 


5.8% 


2.57% 


Wa 


0.86 


1.21 


0.78 


0.29 


0.07 


ris 


0.65% 


0.66% 


0.65% 


0.62% 


0.42% 


7 


142% 


74% 


29% 


13.4% 


4.6% 


''f 


11% 




4.0% 


1.7% 


0.95% 


b? 




5.1% 


5.5% 


1.8% 


0.93% 


bl 


16% 




4.9% 


2.1% 


1.07% 


b^ 




4.6% 


5.3% 


2.7% 


1.18% 


bl 


22% 




4.5% 


1.6% 


0.92% 


bf 




5.9% 


3.5% 


1.7% 


0.92% 


bl 


26% 




4.7% 


2.0% 


1.05% 


b? 




8.9% 


9.6% 


2.7% 


1.20% 



Table 6. Marginalized error-bars for each parameter and different 
probes. For 200 deg-^, except the last column which is for 5000 
deg2. All with Planck and SN-II priors. 



5.11 Marginalized errors 

Table [6] shows the marginalized 1-sigma error-bars of differ- 
ent parameters for each probe. This corresponds to 200 deg'^ 
with Planck and SN-II priors. The last column is for 5000 
deg^. The table reflects in other parameters what we have 
already found for the different FoMs: that the FxB probe 
provides improvements in errors by factors of a few with 
respect to the errors in separate MAG or RSD probes. In 
particular note how the error in bias goes from ~ 5 — 8% to 
~ 1 — 2%. The error in does not change as is dominated 
by priors. This is fine as our goal is to measure the growth 
evolution and not the overall normalization. 



6 SURVEY COMPARISON 

Table [7] shows a comparison of different surveys and probes. 
In the upper section of Table [7] we compare different WL 
probes for the parent F5000 photometric sample. The fol- 
lowing entries correspond to probes for the brighter spec- 
troscopic samples, smaller area surveys and different probes 
either combined over different areas or cross correlated. 

6.1 Shear-Shear, Galaxy-Shear and MAG 

In the first three entries in Table [T] we compare MAG with 
Galaxy-Shear (GS+GG) and Shear-Shear (SS) using the 
same assumptions, fiducial model and priors for a 5000 deg'^ 
survey. Using SS, Ck^k , is the most conservative use of WL 
as it does not requiere any galaxy biasing modeling. Note 
that in GS+GG we also include the galaxy-galaxy auto- 
correlations, i.e. it includes both Cg.^ and Cg-g^, but does 
not include SS, C^iKj , or galaxy cross-correlation, Cg^^g- , for 
i ^ j- Magnification (MAG) includes both auto and cross- 
correlations, i.e. Cg. ,g- , for all i and j. As expected, the MAG 
resuh for F5000 is lower than SS or the GS+GG, this is due 
to galaxy bias. When bias is fixed (entry #2b in the table). 
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MAG is larger despite the relatively poor radial resolution, 
(Jz = 0.05(1 + z). The galaxy-shear cross-correlation alone, 
i.e. Cgk, has lower noise than MAG, i.e. Cgg, as can be seen 
in the covariance of Eq]60| 



(see S2.3I. We consider two different surveys for a BAO 



C0V{Cgg,Cgg) 
COV(Cgk,Cgk) OC CggCkk 



(95) 
(96) 



where Ckk is much smaller than Cgg. This means that with- 
out noise GS performs better than MAG. But systematic 
effects will be quite different for measuring accurate shapes 
versus accurate photometry. It is also likely that in practice 
there are higher densities of galaxies available for MAG than 
for Shear as we move to higher redshifts and this has not 
been taken into account here. 

An advantage of our approach is that we can include all 
observables in the same FM (with appropriate covariances) 
to provide a joint analysis of MAG-I-GS+GG-I-SS includ- 
ing all covariances (label WL-all in the Table). The WL-all 
case (entry #4 in the Table) includes all these weak lensing 
probes, and of course gives the best results. The improve- 
ment of WL-all over GS-I-GG is over a factor of 6: compare 
entries #4 and #3. This is due to the direct measurement of 
the matter amplitude (ie growth) by SS which helps break- 
ing degeneracies with bias in GS+GG and MAG. 



6.2 Bright spectroscopic samples 

As defined in ij4]the B5000 sample is close to spectroscopic 
as it has very good photo-z, ~ 0.0035(1 + z), but a very 
high completeness (of 10000 galaxies/deg^ to zas < 22.5). 
This is motivated by the PAUCam concept (see §1), but we 
also show results for more standard (future) spectroscopic 
flux limited samples. One of them is called B5000spec: it 
has az = and 1000 galaxies/deg^. Otherwise B5000spec is 
identical to B5000. A second sample is called B5000spec21.5 
is flux limited to iAB < 21.5 (instead of iAB < 22.5 in 
B5000spec or B5000) with only 500 galaxies/deg'^. This gives 
a mean redshift of 0.4 and about half the volume of iAB < 
22.5 (see also long dashed line in Fig|4|. As shown in Table 
[t] (comparing entries 6-7), the RSD results for B5000spec 
are lower than for B5000, indicating that density is more 
important than spectroscopic accuracy as found in §5.6| and 
§5.8[ This is a good validation of the PAUCam concept. The 
result for B5000spec21.5 is almost a factor of 2 lower, but 
still quite competitive given that it will be much faster to 
get spectra for galaxies which are 1 magnitude brighter. 

When the B5000 sample is a subset of F5000 (i.e. in 
the FxB case) we can also estimate MAG for B5000 as we 
will have both good redshift accuracy and good photometry 
to do shear or magnification. This is shown as entry #5 in 
Table [7| Note how MAG over B5000 provides slightly better 
FoM„ (but lower FoM^t) than MAG over F5000. When bias 
is fixed MAG over B5000 is always better than MAG over 
F5000 (compare #2b and #5b), despite the larger volume 
and depth in F5000. This of course is due to the better 
resolution in the radial bins. 

BAO is also included in our forecasts for comparison. 
We follow the method of Seo & Eisenstein (2007) using the 
same fiducial values and parameters as for the other cases 



follow-up of F5000. The usual B5000, which for FxB can be 
considered a brighter subset of F5000 and therefore lies in its 
foreground, with most of the redshifts at z < 0.5. the other 
is B5000highz which corresponds to a selection (e.g. emis- 
sion line galaxies) which is biased towards higher redshift 
targets out of F5000, with around 2000 galaxies/deg^, used 
to maximized the BAO probe. We use a constant number of 
target galaxies per unit redshift in the range 0.2 < z < 1.7 
which has about 2.5 larger volume than B5000. The BAO 
result for B5000highz is about a factor of 2 better than for 
B5000, as expected. This is about the largest BAO FoM 
that one could get based on targets from the parent F5000 
sample as we have assumed it covers the same volume as 
F5000 and has quite a large density. The result is almost as 
large as the FoM^, for MAG in B5000 (entry #5), but note 
that with MAG we can also measure 7, which is not possible 
with BAO. As mentioned earlier, doing MAG over B5000 is 
in fact only a subset of what we can do when we combine 
the good radial resolution in B5000 with the photometrical 
(and WL) information in F5000, i.e. the FxB combinations 
which will be discussed next. 



6.3 Smaller area surveys 

In the 11th- 13th rows of Table [7] we see how a modest sur- 
vey, with only 200 deg^ can give a better FoM based on cross- 
correlation (FxB200 in the table) than much larger surveys 
(F5000 or B5000) that consider only either MAG, SS, GS, 
RSD or BAO separately, i.e. all entries above in Table [7j 
except entries #3 and #4 which combine several probes. In 
particular, the BAO PAU Survey proposed in Benitez et al. 
(2009) with 5000 deg^ yields a FoM„ = 47 (entry #9 in the 
Table) which is a factor of 2-3 times smaller than the com- 
bination of RSD with MAG (or WL probes) over 200 deg^ 
(entries #11-13). This represents a factor of 50-70 improve- 
ment in the ratio of benefit over cost, if we assimilate cost 
with area and benefit with FoM. In addition, the RSD+WL- 
all cross-correlation also provides constraints on 7 of about 
5% accuracy, which cannot be obtained from BAO. 

If we combine the Fisher matrix of MAG and RSD from 
two large surveys (F-I-B5000) we do better than FxB200, 
but not by far. When we include all WL and RSD in the 
F-fB5000 combination (row #18-19) the improvement is 
more significant, overcoming the FxB200 result by a fac- 
tor of 3- 7. A s the combined FoM scales roughly with the 
area (Eq87l, doing a survey with ~ 800 deg'^ with the 
FxB cross-correlation will provide similar combined FoM 
than the combination of the two larger but separate sur- 
veys (F5000-I-B5000), which correspond to a total of 10000 
deg^. This represents a cost saving of about a factor of 10, 
if we assume that cost scales with survey area. Also recall 
that F-l-B only uses RSD over B5000 and lensing over F5000, 
while FxB include all probes and all cross-correlations. 

In the case of the PAU Survey, we find from photo-z 
simulations that the narrow-band filters also improve the 
photo-z accuracy of the Faint sample from 0.05(l-|-z) to 
about 0.03(H-z). If we include this improvement, the FoM™^ 
for FxB200 increases by another 30%. 
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Probe 


Sample 




Volume 


crz/(l + 2) 


FoM™ 


FoM-y 


FoM„.y 


# 






Complete 


Gpc^/h^ 


Photo-z error 


7 free 




xlO^ 






Photometric 


















22.5 < iAB < 24 














1 


Shear-Shear (SS) 


F5000 


1/2 


18 


0.05 


84 


39 


3.2 


2 


MAG 


F5000 


1/2 


18 


0.05 


54 


6 


0.33 


2b 


MAG (bias fixed) 


F5000 


1/2 


18 


0.05 


142 


59 


8.4 


3 


Galaxy-Shear (GS-I-GG) 


F5000 


1/2 


18 


0.05 


91 


29 


2.6 


4 


Weak Lensing (WL-all) 


F5000 


1/2 


18 


0.05 


331 


47 


16 






Spectroscopic 


















■Lab < 22.5 














5 


MAG 


B5000 


1/2 


7 


0.0035 


64 


0.9 


0.06 


5b 


MAG (bias fixed) 


B5000 


1/2 


7 


0.0035 


196 


85 


17 


6 


RSD 


B5000 


1/2 


7 


0.0035 


27 


17 


0.47 


7 


RSD 


B5000spec 


1/20 


7 





20 


17 


0.35 


8 


RSD 


B5000spec21.5 


1/26 


3.5 





10 


16 


0.2 


9 


BAO 


B5000 


1/2 


7 


0.0035 


47 






10 


BAO 


B5000highz 


1/10 


18 





78 






Small: 200 deg^ 


11 


MAG-I-RSD 


FxB200 


1/2 


0.7 


0.05-0.0035 


86 


21 


1.8 


lib 


MAG+RSD (bias fixed) 


FxB200 


1/2 


0.7 


0.05-0.0035 


197 


39 


7.7 


11c 


SS-fRSD 


F+B200 


1/2 


0.7 


0.05-0.0035 


37 


18 


0.7 


12 


GS-I-GG-I-RSD 


FxB200 


1/2 


0.7 


0.05-0.0035 


105 


21 


2.2 


13 


WL-all-l-RSD 


FxB200 


1/2 


0.7 


0.05-0.0035 


198 


26 


5.2 


Combine Independent 


14 


BAO-I-RSD 


B5000highz 


1/10 


18 





100 


27 


2.7 


15 


BAO-I- WL-all 


F+B5000highz 


1/2-1/10 


7-18 


0.05-0 


384 


48 


17 


16 


BAO-I-RSD-I- WL-all 


F+B5000highz 


1/2-1/10 


7-18 


0.05-0 


597 


66 


40 


17 


MAG+RSD 


F+B5000 


1/2 


7-18 


0.05-0.0035 


82 


29 


2.4 


17b 


MAG+RSD (bias fixed) 


F+B5000 


1/2 


7-18 


0.05-0.0035 


686 


112 


77 


17c 


SS+RSD 


F+B5000 


1/2 


7-18 


0.05-0.0035 


262 


58 


15 


18 


WL-all+RSD 


F+B5000 


1/2 


7-18 


0.05-0.0035 


571 


66 


38 


19 


SS-RSD 


F+B5000spcc 


1/2-1/20 


7-18 


0.05-0 


257 


58 


15 


Cross-correlation 


20 


MAG+RSD 


FxB5000 


1/2 


7-18 


0.05-0.0035 


645 


46 


30 


20b 


MAG+RSD (bias fixed) 


FxB5000 


1/2 


7-18 


0.05-0.0035 


1141 


139 


159 


20c 


SS+RSD 


FxB5000 


1/2 


7-18 


0.05-0.0035 


316 


60 


19 


21 


WL-all+RSD 


FxB5000 


1/2 


7-18 


0.05-0.0035 


2879 


87 


251 


22 


WL-all+RSD 


FxB5000spec 


1/2-1/20 


7-18 


0.05-0 


2113 


74 


159 


23 


WL-all+RSD 


FxB5000spec21.5 


1/2-1/26 


3.5-18 


0.05-0 


1509 


65 


98 



Table 7. FoM comparison for different surveys and probes: Shear-Shear (SS), Galaxy-Shear (GS+GG, which includes Galaxy-Galaxy 
autocorrelations), Magnification (MAG), redshift space distortions (RSD) and Baryon Acoustic Oscillations (BAO). The WL-all case 
includes MAG, SS, GG and GS. All cases are marginalized over bias (except for #2b, 5b, lib, 17b and 20b, where bias is fixed) and 6 
cosmological parameters (see Eq|9|. All samples are 5000dcg2 expect for #11-13 which are 200deg2. The 5th column gives the comoving 
volume to the highest redshift used in each analysis. BAO alone (#9-10) does not provide FoM-,. or FoMu,^. 



6.4 Spectroscopic and Photometric combinations 

Entries #20-21 of Table [T] shows that the FxB5000 combi- 
nation are about 10 times larger FoM„^ than F+B (#17- 
18), even when F+B correspond to twice the area (as they 
come from independent surveys and therefore have twice 
as many independent modes). Part of this gain comes from 
the fact that in the FxB combination we have two probes 
(WL and RSD) to measure the bias of the B population, 
while for the F+B combination the bias in B is only measure 
with RSD. The improvement on going from F+B to FxB is 
larger in FoM^ than in FoM.^. Note that the gain is smaller 
for SS+RSD combination (compare #17c to #20c) because 
bias is already well measured with F+B with Shear. This 
can also be seem comapring the cases with fixed bias (#17b 
and #20b) . When bias is fixed in the FxB case (#20b) 



there is still a significant gain with respect to the free bias 
case (#20), indicating that there is still room for important 
benefits if we can measure bias in some other ways {[ 2.41. 



In principle one can also infer the bias of the B spectro- 
scopic population by identifying the same bright population 
within the photo-z {iAB < 24) sample and doing the cross- 
correlation of the faint and bright galaxies. This can improve 
the F+B result by a factor of 5, but is still a factor of 2 lower 
than FxB because of the larger photo-z errors. 

This gain is illustrated in Fig |13[ which shows the 3 pro- 
jections {wojWa,"/) in the FoM^^, ellipsoid. Alone, the RSD 
results (purple) in the spectroscopic sample (BSOOOspec) 
are less powerful in constraining FoMuj^ than Shear-Shear 
(blue), but the improvement resides in the complementar- 
ity of the combination to break degeneracies. This can be 
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Figure 13. Contours (Ax^ = 2.3) of wq — Wa — 7 in weak lensing (only Shear-Shear) over a Faint sample with iab < 24 (F5000, blue) 
as compared with RSD over a bright spectroscopic subset (BSOOOspec, purple) targeting 1000 galaxies per square degree to i^^g < 22.5. 
Bias evolution {b{z)) and other cosmological parameters are marginalized over. The red contours correspond to the combination of both 
probes over different sky areas (F-l-B), while the yellow contours show the combination over the same area (FxB), including also all 
cross-correlations between shear and galaxies in both samples. This produces over a factor 10 larger combined FoM. 



seen in the red contours (F-|-B5000spec) which combine both 
RSD and Shear-Shear assuming they are independent. When 
the BSOOOspec sample is a subset over the same area of the 
F5000 sample, we can also include all their cross-correlations 
and provide a second route to measure bias evolution. This 
benefits both WL and RSD, and the combination is over a 
factor of 10 larger (yellow FxB contours in the figure) than 
F+B and a factor of 10 times larger than all WL in the Faint 
sample alone (entry #4 in Table [7|. 

Figure |5] shows the corresponding result using MAG 
and RSD over 200 deg^, instead of Shear and RSD over 5000 
deg^. In this case the gain in the combination FxB (see Table 
[3| is over a factor of 10 times larger than F-l-B and over a 
factor 50 for the one in MAG in the faint sample alone. The 
extra benefit here comes from the higher sampling rate in 
the PAUCam strategy. 

The BAO-f RSD (B5000highz) combination (entry #14) 
is comparable to the MAG-f RSD in F-f B5000 (#17) or the 
FxB200 combinations (#11). BAO and RSD have also been 
used to constrain bias (e.g. see Amendola, Quercellini and 
Giallongo 2005, Cinzia, Amendola and Branchini 2011). But 
this is a factor of 5 lower than the WL-all-f RSD independent 
combinations (#18) and a factor of 90 lower than the FxB 
cross-correlation result (#21). The results for B5000spec21.5 
(last entry in the Table) are about a factor of 1.6 lower in 
the combined FoM than the corresponding iab < 22.5 result 
(#22). 



6.5 Spectroscopic follow-up strategy 

The WLxRSD (or MAGxRSD) strategy appears much bet- 
ter than following a BAO approach over the same B5000 



sample. The combination of WL and RSD produces a FoM^, 
about 24 times larger that BAO alone over B5000. The rea- 
son is that in BAO we are ignoring the P{k) amplitude 
and WL distance ratios information. In fact, RSD alone can 
give comparable FoM than BAO despite the additional bi- 
asing parameters. It also produces interesting constraints in 
7 which BAO cannot measure. 

A potential advantage of a BAO follow-up would be to 
sample larger volumes as in B5000highz which corresponds 
to a selection of higher redshift targets out of F5000, with 
around 2000 galaxies/deg^ with constant density per unit 
redshift between 0.2 and 1.7. This is shown in the #10 entry 
of Table [7| The improvement in FoM^, is significant but still 
far from WL alone (#4). The combination of BAO and RSD 
produces only a modest improvement over the BAO alone 
result, similar to MAG-I-RSD from separate surveys. 

To have a fair comparison or optimization of which 
could be the best strategy (going deep as in B5000highz 
or shallower as B5000spec) we would need to add BAO to 
WL and RSD as separate probes (F-|-B). This is not to- 
tally straight forward as these probes could be correlated. 



Our WL forecast already includes BAO wiggles (see S5.9I. 
What we can do here is to combine the BAO measurement 
with WL and RSD, as if they were from independent sam- 
ples. This is partially true as there is little volume overlap 
between B5000highz and B5000spec (where the lenses to 
F5000 reside). In any case, this is an upper bound of the 
true result, because the covariance between WL and BAO 
will lower the FoM. Results are shown as entries #15-16 in 
Table [7] The BAO-l-WL-all combination (entry #15) is sim- 
ilar to SS-fRSD (entry #19). When we combine everything 
BAO-l-WL-alH-RSD (entry #16) the outcome is of course 
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higher. But it is still over a factor of 6 lower than the FxB 
WL-all+RSD cross-correlation {^21). This indicates that a 
shallower spectroscopic follow-up strategy provides better 
returns in terms of FoM. The legacy value of such spectro- 
scopic follow-up could also be larger as we will be able to 
relate dark matter, as traced by WL, to galaxy formation, 
as traced by galaxies over the same structures. 



7 CONCLUSION 

In this paper we advocate the use of the cross-correlation of a 
deep (faint) photometric sample with a foreground brighter 
spectroscopic subset using narrow radial redshift bins. The 
foreground sample could also be photometric with very good 
photo-z of 0-^/(1 + 2) ~ 0.0035. We consider three differ- 
ent types of probes in the analysis: 1) angular clustering 
from galaxy-galaxy autocorrelation in narrow redshift bins, 
2) weak lensing (WL) from shear-shear, galaxy-shear and 
magnification (i.e. galaxy-galaxy cross-correlation), 3) red- 
shift space distortions (RSD), from the ratio of transverse to 
radial modes. The combination of such measurements pro- 
vides a significant improvement in the forecast for the evo- 
lution of the dark energy equation of state, given by w{z), 
and the cosmic growth evolution, given by 7. 

This improvement comes from the use of narrow red- 
shift bins and the measurement of galaxy bias, which affects 
both RSD and WL cross-correlations, but in different ways. 
For transverse modes, biasing can be obtained by compar- 
ing WL clustering with galaxy angular clustering, but this 
is also degenerate with geometrical factors that depend on 
cosmic expansion evolution. In redshift space, bias can be 
measured directly by comparing line-of-sight to transverse 
clustering. The combination of both methods results in a 
very accurate determination of bias evolution, also breaking 
the degeneracy of growth with cosmic expansion evolution. 
We have argued that this gain is only fully accomplished 
when both WL and RSD measurements are done over the 
same area of the sky, which allows better measurement of 
bias and the full understanding of cross-correlations. Once 
bias is known, the foreground sample gives much better con- 
straints than the deeper photometric sample, thanks to its 
better radial resolution. 

A new figure of merit, FoM^j, that encodes both cos- 
mic expansion 'w(z) and growth evolution 7 is introduced to 
explore changes with the sample area. A, the redshift com- 
pleteness, -q, the magnitude limit of the bright sample, mi, 
the photo-z or spectroscopic error, and the uncertainty 
in the fraction of photo-z transitions Ar defined in §3.12| 
For the magnification (MAG) and RSD combination over 
the same area (FxB) we find 



we have explored, i.e. A £ [1, 50], > 0.01, m; G [21.5, 23.0], 
(Jz < 0.7 and < 0.6. If we include all WL probes: shear- 
shear, galaxy-shear and magnification, the above FoM^-y is 



FoM^T- ~ 2080 A°-^^ jf-'^ 1.26 



mi —22.5 



where A is the survey area in units of 200 deg^, m; is the 
magnitude limit of the bright sample, is the photo-z error 
of the bright sample in units of 0.01(1 -I- 2) and is the 
absolute uncertainty in the fraction of photo-z transitions in 
units of 0.18. This equation can be used to extrapolate the 
predictions to different survey parameters within the ranges 



about a factor of 2.9 (8.6) times larger than Eq 97 for 200 
(5000) deg^ These results are robust and conservative. We 
only use linear scales (fcmin = 0.1) and a relatively large red- 
shift bin-width of A2 ~ 0.014(1 -I- 2). Better results can in 
principle be obtained using smaller scales and narrower red- 
shift bins, but this requieres further assumptions and mod- 
eling. We use a conservative modeling of biasing, which uses 
four parameters for each population (i.e. faint and bright 
subsamples) and no priors. With two bias parameters per 
population the above FoM increases by about 40%. If we fix 
bias the improvement is 4.3 (5.3) for 200 (5000) deg^. 

Shear-Shear alone (excluding galaxy correlations) pro- 
duces smaller FoM than magnification alone with bias fixed, 
but 10 times lager when bias is not known. Galaxy-Shear 
(including galaxy-galaxy autocorrelation and galaxy bias) 
can produce a factor 8 times better FoM^ than MAG 
(i.e. galaxy-galaxy cross and auto-correlation) because the 
latter is subject to larger sampling error from the back- 
ground galaxy-galaxy autocorrelation. All WL probes ben- 
efit equally from the boost in the cross-correlation with 
RSD. They provide similar cosmological information (see 
Van Waerbeke 2009), so that the final FoM will depend on 
the size of systematic errors and the number density and 
volume that we can trace in each case. As illustrated by our 
results in Eq]97[ magnification alone, without shape infor- 
mation, can be used to provide very competitive cosmolog- 
ical information. 

We have not considered here photo-z bias or shear mul- 
tiplicative or additive biases as they are more survey de- 
pendent and have been studied elsewhere (see Huterer et al. 
2006, Bernstein 2009 and references therein). But we have 
found new requirements on the uncertainties of photo-z tran- 
sitions for the magnification, galaxy-shear and shear-shear 
probes to work. For all cases, the transitions Tij (or photo-z 
contamination) needs to be known to about 1% (absolute er- 
ror in Tij) if we want only a small degradation in the FoM™^ 



(see Fig 12 and Eq 94 1. These values look within reach for 
future surveys. But note that as indicated by Eq |93| this 
requirement increases as (Hax where Imax is the maximum 
multipole (or smaller angular scale) used in the fit. In our 
case we only use linear scales (Eq |13| |. A fit using non-linear 
scales (Imax — 6000) needs to know photo-z transitions 
to 0.4% and also requieres understanding of how baryonic 
physics could change the non-linear P{k) (Van Daalen et al. 
2011, Semboloni et al. 2011). 

Given current limitations on systematics of shear mea- 
surements, specially for deeper ground based surveys, push- 
ing for magnification as a cosmological tool is a good alter- 
native to implement the WL-I-RSD cross-correlation advo- 
cated here. A clear advantage for magnification would be to 
consider deeper photometric samples (such as LSST which 
can reach r ~ 27.5, Ivezic, et al. 2011), where galaxy shapes 
are dominated by atmospheric noise but we can still measure 
magnitude and photo-z errors (see Van Waerbeke et al. 2010, 
Namikawa, Okamura & Taruya 2011). In this case, the mag- 
nification and RSD cross-correlation presented in this paper 



© 2011 RAS, MNRAS 000, 1-10 



Cosmology from lensing & redshift distortions 29 



could play a significant role, provided we can perform an 
appropriate deep spectroscopic follow-up of the foreground 
lenses and provided we can calibrate the fraction of photo-z 
transitions. The PAUCam concept, based on doing photo- 
metric redshifts from 40 narrow (lOOA) and 6 broad bands, 
seems well suited for such deep spectroscopic follow-up as 
we have found that the characteristic PAU photo-z error of 
= 0.0035(1 + z) only degrades the FoM„^ by ~ 10% (e.g. 
FigP ^ 



11 



and Eq 97 1 



To design a spectroscopic follow-up survey of a pho- 
tometric (WL) survey, we find that we can obtain 6 times 
better FoM^-y by targeting the foreground lenses to do cross- 
correlations, as proposed here, than trying to go deeper to 
measure BAO (compare entry #21 to #16 in Table 6), as has 
been proposed in several future surveys, such as Euclid, Big- 
BOSS, HETDEX or FOCAS (WFMOS). Measuring BAO 
with a 5000 deg^ PAU Survey (Benitez et al. 2009) provides 
2-4 times smaller FoM„ than doing RSD-I-MAG (or WL) 
cross-correlation with the same PAUcam instrument over a 
much smaller 200 deg'^ area (compare #11-13 to #9 in Table 
[7|. This represents a factor 50 to 100 improvement in profit 
over cost, defined as FoM achieved per unit area. 

In our analysis we have only considered the gains in the 
photometric and spectroscopic cross-correlation using the 
2-point clustering statistics. Higher order correlations will 
certainly add additional benefits both in the weak-lensing 
probes (e.g. see Gaztanaga & Bernardeau 1998), in clus- 
tering and in galaxy biasing (Frieman & Gaztanaga 1994, 
Gaztanaga & Scoccimarro 2004). Besides the improvement 
in cosmological inference, the cross-correlations of spectro- 
scopic and photometric surveys should also provide key in- 
formation on the relation between galaxy evolution and dark 
matter growth, which will bring some new light into the 
problem of galaxy formation. 
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